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Abstract Over the last several years, we have embarked on a long term effort to exploit the strong potential that hot B subdwarf 
(sdB) pulsators have to offer in terms of asteroseismology. This effort is multifaceted as it involves, on the observational front, 
the acquisition of high sensitivity photometric data supplemented by accurate spectroscopic measurements, and, on the theoret¬ 
ical and modeling fronts, the development of appropriate numerical tools dedicated to the asteroseismological interpretation of 
the seismic observations. In this paper, we report on the observations and thorough analysis of the rapidly pulsating sdB star (or 
EC14026 star) PG 1219-t534. Our model atmosphere analysis of the time averaged optical spectrum of PG 1219+534 obtained 
at the new Multiple Mirror Telescope (MMT) leads to estimates of Teff = 33,600 ± 370 K and logg = 5.810 ± 0.046 (with 
log V(He) /N(H) = -1.49 ± 0.08), in good agreement with previous spectroscopic measurements of its atmospheric parameters. 
This places PG 1219+534 right in the middle of the EC14026 instability region in the logg - Tcti plane. A standard Fourier 
analysis of our high signal-to-noise ratio Canada-France-Hawaii Telescope (CFHT) light curves reveals the presence of nine 
distinct harmonic oscillations with periods in the range 122 - 172 s, a significant improvement over the original detection of 
only four periods by Koen et al. (1999). On this basis, we have carried out a detailed asteroseismic analysis of PG 1219+534 
using the well-known forward method and assuming that the observed modes have £ < 3. Our analysis leads objectively to 
the identification of the (k, £) indices of the nine periods observed in the star PG 1219+534, and to the determination of its 
structural parameters. The periods all correspond to low-order acoustic modes with adjacent values of k and with C = 0, 1, 2, 
and 3. They define a band of unstable modes, in close agreement with nonadiabatic pulsation theory. Furthermore, the average 
dispersion between the nine observed periods and the periods of the corresponding nine theoretical modes of the optimal model 
is only ~ 0.6%, comparable to the results of a similar analysis carried out by Brassard et al. (2001) on the rapid sdB pulsator 
PG 0014+067. On the basis of our combined spectroscopic and asteroseismic analysis, the inferred global structural parameters 
of PG 1219+534 are Gff = 33,600 ± 370 K, logg = 5.8071 ± 0.0057, logMenv/M. = -4.254 ± 0.147, M. = 0.457 ± 0.012 Mq, 
R/Ro = 0.1397 ± 0.0028, and L/Lq = 22.01 ± 1.85. Combined with detailed model atmosphere calculations, we estimate, in 
addition, that this star has an absolute visual magnitude My = 4.62 ± 0.06 and is located at a distance rf = 531 ± 23 pc (using 
V = 13.24 ± 0.20). Finally, if we interpret the absence of fine structure (frequency multiplets) as indicative of a slow rotation 
rate of that star, we further find that PG 1219+534 rotates with a period longer than 3.4 days, and has a maximum rotational 
broadening velocity of V sini <2.1 km.s^'. 
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1. Introduction 

Hot subdwarf B (sdB) stars are commonly interpreted as low- 
mass (~ 0.5 Mo) core helium burning objects that belong to 
the so-called Extreme Horizontal Branch (EHB). Having ex¬ 
tremely thin and mostly inert hydrogen rich residual envelopes 
(Menv < 0.02 Mq), these stars remain hot (Eeff ~ 22,000 - 
42,000 K) and compact (logg ~ 5.2 - 6.2) throughout their en¬ 
tire EHB lifetime and eventually evolve toward the white dwarf 
cooling sequence without experiencing the Asymptotic Giant 
Branch (AGB) and Planetary Nebulae (PN) phases of stellar 
evolution (see Dorman et al. 1993). If the fate of sdB stars is 
relatively well understood, their origin - i.e., how they form - 
is still unclear, however. The main difficulty resides in identify¬ 
ing how sdB stars manage to remove all but a very small frac¬ 
tion of their envelope at exactly the same time that they reach 
a sufficient core mass for helium flash ignition at the tip of the 
Red Giant Branch (RGB). Various formation scenarios have 
been proposed over the years to solve this problem, involving 
either single star evolution with enhanced mass loss at the tip 
of the RGB (D’Cruz et al. 1996), or binary evolution including 
various concurrent channels from common envelope ejection, 
stable Roche lobe overflow, to the merger of two helium white 
dwarfs (see, e.g., Han et al. 2002, 2003 and references therein). 
Interestingly, these authors suggest that the latter binary evolu¬ 
tion scenarios would lead to a somewhat broader dispersion of 
stellar masses among sdB stars (between 0.30 - 0.70 Mq) than 
is currently assumed from canonical EHB models. 

Renewed interest in this phase of stellar evolution has come 
in part from the recent discoveries of two new classes of mul- 
tiperiodic nonradial pulsators among sdB stars which enable, 
in principle, the asteroseismic probe of their internal struc¬ 
ture. The rapid sdB pulsators known as the V361 Hya stars, 
but more commonly referred to as the EC14026 stars, were 
the first to be observationally detected by colleagues from the 
South African Astronomical Observatory (SAAO; Kilkenny 
et al. 1997). Interestingly, at the same epoch, their existence 
was also independently predicted from theoretical considera¬ 
tions based on the identification of an efficient driving mecha¬ 
nism for the pulsations (Charpinet et al. 1996). The mechanism 
uncovered is a /r-effect due to local overabundances of heavy 
metals - especially of iron - produced by microscopic diffu¬ 
sion processes which are at work in the envelope of these stars 
(Charpinet et al. 1997). Various surveys by the SAAO group, by 
the Montreal group (see Billeres et al. 2002), by a Norwegian- 
German-Italian team (see, e.g., Silvotti et al. 2002), and by 
other teams (see the reviews of Charpinet 2001 and Kilkenny 
2002) rapidly led to additional discoveries, bringing the total of 
known EC14026 pulsators to 33 at the time of this writing. The 
EC 14026 stars tend to cluster at high effective temperatures 
and surface gravities (near ~ 33,500 K and logg ~ 5.8), 
although outliers indicate a relatively large dispersion of these 
stars in the log g - plane. The periods detected in EC 14026 
stars typically range from ~ 100 s to ~ 200 s, but can be sub¬ 
ported here were obtained at the MMT Observatory, a joint facility of 
the University of Arizona and the Smithsonian Institution. 

** This study made extensive use of the computing facilities offered 
by the Calcul en Midi-Pyrenees (Calmip) project, France. 


stantially longer for class members having lower surface gravi¬ 
ties (e.g., 290-600 s for PG 1605+072; see Charpinet 2001 and 
references therein). Their mode amplitudes usually span a rel¬ 
atively wide range, but typical values are of several millimag- 
nitudes. In most cases, the short periods are entirely consistent 
with low-order, low-degree radial and nonradial acoustic waves 
(Stobie et al. 1997; Charpinet et al. 1997, 2001). However, in 
low surface gravity sdB pulsators, the presence in the observed 
period spectrum of low-order g-modes or mixed modes having 
both p- and g-mode properties cannot be excluded and some¬ 
what complicates the identification of the true nature of the 
modes being detected (Charpinet et al. 2002). Remarkable sim¬ 
ilarities exist between appropriate EHB pulsation models (i.e., 
our so-called “second generation” models which take radiative 
levitation of iron into account) and the most basic properties of 
EC14026 pulsators (see the review of Charpinet et al. 2001). 

The long period sdB variables (the PG1716+426 stars, but 
sometimes also referred to as the “Betsy stars” or the “IpsdBV 
stars”) were discovered only a couple of years ago (Green et al. 
2003). PG1716+426 stars oscillate on much longer timescales 
than their EC14026 counterparts, with periods ranging from 
~ 45 minutes to ~ 2 hours. This implies that relatively high 
order g-modes are involved. Moreover, unlike the EC 14026 
pulsators, these stars populate the low-temperature/low-gravity 
domain of the logg - Teff plane where sdB stars are found. 
Remarkably, however, the same mechanism responsible for the 
oscillations in the EC 14026 stars is thought to operate in the 
long period sdB pulsators as well, but destabilizing this time 
high-order, f > 3 gravity modes (Eontaine et al. 2003). 

Eor the asteroseismologist, the existence of these two dis¬ 
tinct classes of pulsators among the hot B subdwarfs is of par¬ 
ticular interest, as p-modes and g-modes probe different re¬ 
gions of their interiors (see Charpinet et al. 2000). One can 
therefore hope to gather complementary information on EHB 
stars in general from the study of these two classes of pul¬ 
sators. While it is too early to assess the real asteroseismolog- 
ical potential of the newly discovered PG1716+426 stars, the 
EC 14026 pulsators have already proved to be excellent labo¬ 
ratories for asteroseismic studies. In the exploratory work of 
Brassard et al. (2001), the acquisition of high signal-to-noise 
(S/N) ratio data at the Canada-Erance-Hawaii 3.6 m Telescope 
(CEHT) combined with efforts to develop a new objective 
global optimization method for asteroseismology led to the first 
successful attempt to match all the detected periods of a short 
period sdB pulsator. These authors succeeded in reproducing 
simultaneously (at the ~ 1% level) 13 periods extracted from 
the light curve of the star PG 0014+067. This led to the first 
asteroseismological determination of the fundamental param¬ 
eters of an sdB star, such as its effective temperature, its sur¬ 
face gravity (to a greatly improved accuracy), its total mass, 
and the mass of its H-rich envelope. The derivation of that lat¬ 
ter quantity, which is of utmost importance in the context of 
Extended Horizontal Branch (EHB) stellar evolution to which 
sdB stars are bound, is a pure product of asteroseismology and 
currently cannot be measured by any other known technique. 
Brassard et al. (2001) have emphasized the important fact that 
their asteroseismological solution is entirely consistent with 
both the atmospheric parameters (Teff and logg) derived from 



Charpinet et al.: Structural Parameters of the Hot Pulsating B Subdwarf PG 1219+534 from Asteroseismology 


3 


independent quantitative spectroscopy and the mode stability 
predictions obtained from the nonadiabatic pulsation calcula¬ 
tions based on the so-called second generation models. Such a 
consistency between three independent aspects of the modeling 
of these stars - namely model atmosphere computations, pe¬ 
riod distribution calculations, and nonadiabatic mode stability 
predictions - certainly adds to the credibility of that asteroseis- 
mological result. 

Since the work of Brassard et al. (2001), we have greatly 
perfected our technique and successfully extended such studies 
to several pulsating hot B subdwarfs, as well as several white 
dwarf pulsators (see, e.g., Charpinet et al. 2003; Brassard & 
Fontaine 2003). In this paper, we report on the case of the hot 
B subdwarf PG 1219+534 (aka V* KY UMa), a star first rec¬ 
ognized as an sdB pulsator by Koen et al. (1999). These au¬ 
thors have identified four periodicities in the discovery light 
curve of this star. In this study, we more than double that num¬ 
ber by reaching a total of nine independent periods detected, 
based on high S/N ratio light curves obtained during dedicated 
observations at the CFHT. This significant improvement at the 
level of the photometry added to new improved spectroscopy 
on that star allows us to perform a full asteroseismic analysis 
of PG 1219+534. This ultimately leads to the determination of 
its fundamental parameters, including its total mass and H-rich 
envelope mass, i.e., two key parameters for constraining for¬ 
mation and evolution scenarios of sdB stars in the future. 

In the following sections, we fully report on our combined 
photometric and spectroscopic observations of PG 1219+534 
(Section 2) and discuss the detailed analysis of the photometric 
light curves obtained for that star (Section 3). Then, we present 
and discuss in detail the outcome of the thorough asteroseismic 
analysis we have carried out on PG 1219+534 based on these 
new observations (Section 4). A summary and conclusions up¬ 
rising from this analysis are then provided in Section 5. 

2. Observations 

2.1. Medium Resolution Spectroscopy 

A reference to PG 1219+534 appears for the first time in the 
Palomar-Green catalog of UV excess objects (Green et al. 
1986). Loosely classified at that time as a subdwarf star with 
no indication of its subclass, the available Johnson photometry 
from Green et al. (1986; B = 12.41) and Stromgren photome¬ 
try from Wesemael et al. (1992; y = 13.235, (b -y) — -0.157, 
(u - b) - -0.265, ml = 0.097) indicate however a relatively 
bright star with colors typical of an isolated sdB object. 

An estimate for the surface parameters of PG 1219+534 
based on the existing absolute multicolor photometry from 
the literature and on dedicated spectroscopy from the 4.2 m 
William Herschel Telescope (WHT) with the spectrograph ISIS 
was first proposed in Koen et al. (1999). Fitting the Balmer 
lines of the spectrum with profiles calculated from a grid of 
LTE H/He model atmospheres (Saffer et al. 1994), these au¬ 
thors adopted the values - 32,800 ± 300 K and logg = 
5.76 ± 0.04. However, Koen et al. (1999) could not provide a 
direct fit of the helium lines and relied instead on comparisons 
with spectra of other sdB stars to estimate a surface helium 


abundance of log A(He)/A(H) ~ -1.18, which is among the 
highest ratios encountered for sdB stars. This rough estimate 
adds additional uncertainty on the T^g and logg values thus 
derived. 

Heber et al. (2000) provided a detailed analysis of a Keck 
HIRES spectrum of PG 1219+534 using line blanketed NLTE 
and LTE model atmospheres to derive the atmospheric param¬ 
eters, metal abundances, and rotational velocity of that star. 
These authors used several indicators based on either hydro¬ 
gen Balmer lines and neutral He lines fitting or He ionization 
equilibrium, but could not reconcile all measurements. In par¬ 
ticular, they pointed out to a “helium problem”, i.e., the inabil¬ 
ity to fit simultaneously well the hydrogen Balmer lines, the 
neutral helium lines, and the Hell 4686 line seen in the spec¬ 
trum of PG 1219+534. This problem seems widespread among 
the hotter sdB stars. By forcing a better fit to the helium lines, 
in particular to the Hell 4686 line, Heber et al. (2000) found 
that the solution would now suggest higher effective tempera¬ 
tures and higher surface gravities than those inferred in their 
global fits (see their Table 5 and Eigure 5). In order to reflect 
these uncertainties, the authors finally adopted values of the at¬ 
mospheric parameters for PG 1219+534 with relatively large 
error bars, namely Teff =34,300^j qqq K, logg = 5.95 ±0.1, and 
log A(He)/A(H) = -1.50 + 0.10. The very sharp metal lines 
seen in the Keck HIRES spectrum also allowed these authors 
to constrain the projected rotational velocity of PG 1219+534 
by setting an interesting limit of V sin / < 10 km s“', indicat¬ 
ing that this star is either a slow rotator, an object seen nearly 
pole-on, or both. 

We must point out here that, at the outset, asteroseismol¬ 
ogy rules out a surface gravity as high as log g = 5.95 for PG 
1219+534 for the simple reason that there are not enough the¬ 
oretical modes available in the range of observed periods. Only 
models with lower values of log g have dense enough period 
spectra, a very robust result that reflects the high sensitivity of 
the pulsation periods on that parameter in sdB stars. A clear 
demonstration of this phenomenon is provided in Eigure 3 of 
Eontaine et al. (1998; see also Eigure 5.13 of Charpinet 1999). 
Accordingly, it will not be surprising to find out below that 
models with log g = 5.95 provide extremely poor fits to the pe¬ 
riod data and that things get worse when increasing further the 
surface gravity. Hence, among the possible atmospheric solu¬ 
tions proposed by Heber et al. (2000), only those with the lower 
surface gravities appear compatible with the period data for PG 
1219+534. 

In order to double check on the spectroscopy of the star PG 
1219+534, we obtained additional measurements from which 
we derived independent estimates for its atmospheric parame¬ 
ters. A medium-resolution (~ lA), high S/N ratio (~ 175 per 
pixel, with 3.1 pixels per resolution element) spectrum of that 
star was obtained recently with the blue spectrograph at the new 
6.5 m Multiple Mirror Telescope (MMT). That spectrum covers 
the range from ~ 4000A to ~ 4950A. Another low-resolution 
(~ 8.7A), high S/N ratio (~ 231 per pixel, with 2.6 pixels per 
resolution element) spectrum from Steward’s 2.3 m Telescope 
was also available. The latter covers the range from ~ 3615A 
to ~ 6900A. Such observations are part of a global program 
aimed at improving the spectroscopic characterization of sdB 
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Figure 1. Model fit (solid curve) to the hydrogen Balmer lines and 
helium lines available in our time averaged high S/N ratio, medium- 
resolution optical spectrum of PG 1219+534 

stars, and in particular of EC14026 stars and long period sdB 
pulsators (the PG 1716+426 stars). 

An integral part of that program is the development of a 
bank of model atmospheres and synthetic spectra suitable for 
the analysis of the spectroscopic data. To this end, we have so 
far computed two detailed grids (one in LTE and the other in 
NLTE) with the help of the public codes Tlusty and Synspec 
(Hubeny & Lanz 1995, Lanz & Hubeny 1995). Each grid is de¬ 
fined in terms of 11 values of the effective temperature (from 
20,000 K to 40,000 K in steps of 2,000 K), 10 values of the 
surface gravity (from logg of 4.6 to 6.4, in steps of 0.2 dex), 
and 9 values of the helium-to-hydrogen number ratio (from 
log A(He)/A(H) of -4.0 to 0.0, in steps of 0.5 dex). These 
grids were originally developped to analyze our MMT data 
and, therefore, our current synthetic spectra are limited to the 
range from 3500A to 5800A, but this can easily be widened 
as needed. We are planning to include metals in the near fu¬ 
ture. More details about these models will be found in Green, 
Eontaine, & Chayer (in preparation). 

We used both grids of synthetic optical spectra to analyze 
our two spectra of PG 1219+534. Using a procedure very sim¬ 
ilar to that prescribed by Saffer et al. (1994), we performed a 
simultaneous fit of the available Balmer lines and helium lines 
within a given grid of profiles from our recent model atmo¬ 
sphere calculations. Eor the NLTE grid, this procedure led to 
the solution for the MMT spectrum shown in Eigure 1, with 
errors given on the parameters corresponding to formal errors 
of the fit. More realistically, based on spectra obtained on dif¬ 
ferent nights to estimate external errors, we adopted - 



Figure 2. Model fit (solid curve) to the hydrogen Balmer lines and 
neutral helium lines available in our time averaged high S/N ratio, 
low-resolution optical spectrum of PG 1219+534 

33,600 + 370 K, logg = 5.810 + 0.046, and logA(He)/A(H) = 
-1.49 + 0.08. Figure 1 indicates that we get a quite decent fit 
to all of the available Balmer and helium lines, except for the 
Hell 4686 line which shows a flux deficit compared to our best- 
fit model. This confirms the finding of Heber et al. (2000) and 
points to a shortcoming in current models. Maybe the inclusion 
of metals in future models could cure this problem. In any case, 
since ours is a global fit, and since the relatively weak Hell 
4686 line plays only a minor role in the overall picture, we like 
to think that ours are realistic estimates of the atmospheric pa¬ 
rameters of PG 1219+534. We note in this context that our LTE 
solution leads to a fit of comparable quality (not shown here) 
but with slightly different parameters: Teff = 33,120 ± 110 K, 
logg = 5.80 ± 0.02, and log A(He)/A(H) = -1.49 + 0.02 (for¬ 
mal errors). This indicates that NLTE effects on the populations 
of hydrogen and helium ions are small in the atmosphere of PG 
1219+534. Of course, of the two solutions, the NLTE one is to 
be preferred because the physics is more accurately treated in 
that case. 

Figure 2 shows our NLTE solution for the low-resolution 
spectrum, with the quoted uncertainties there again referring to 
the formal errors of the fit. More realistic estimates are at least 
double those, and we may adopt Teff = 34040 + 220 K, log g = 
5.86 ± 0.06, and logA(He)/A(H) = -1.51 ± 0.06. Within the 
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Table 1. CFHT High Time Resolution Photometric Observations of PG 1219+534. 


Night 

# 

Run 

Number 

Date 

(UT) 

Start Time 
(UT) 

Start Time 

BID 2453073.5 + 

Sampling 
Time (s) 

Total Number of 

Data Points 

Run length 
(hours) 

1 

cfh-104 

2004-03-09 

06:45 

0.292811 

10 

2221 

6.17 

2 

cfh-105a 

2004-03-10 

06:19 

1.269846 

10 

1336 

3.71 

2 

cfh-105b 

2004-03-10 

10:02 

1.430144 

10 

1471 

4.09 

3 

cfh-106a 

2004-03-11 

07:21 

2.312399 

10 

1610 

4.47 

3 

cfhl06b 

2004-03-11 

11:49 

2.500243 

10 

553 

1.54 

4 

cfh-107 

2004-03-12 

06:14 

3.266342 

10 

3369 

9.36 


quoted uncertainties, both independent measurements are in 
excellent agreement. This is very encouraging and alleviates 
some of the fear (see, e.g., Saffer et al. 1994) that the use of 
only three Balmer lines (in the more limited spectral range of 
our MMT observations) could introduce some systematic ef¬ 
fects. On the basis of a substantial sample of sdB stars observed 
both with the MMT and with the 2.3 m Steward telescope (the 
spectra used here are part of that program), we will establish 
elsewhere that there are indeed no significant systematic differ¬ 
ences between the two approaches. Another extremely encour¬ 
aging development is that, using completely independent tools, 
the referee of this article. Dr. Uli Heber, has kindly commu¬ 
nicated to us that his fit to our low-resolution spectrum using 
comparable models to ours (i.e., NLTE H/He models with no 
metals) leads to T^s - 34200 ± 70 K, logg = 5.86 + 0.01, and 
log A(He)/A(H) = -1.45 + 0.02 (formal errors), in excellent 
agreement with our results. This gives added confidence that 
our estimates of the atmospheric parameters of PG 1219+534 
are secure. 

In practice, in what follows, we will adopt our NLTE so¬ 
lution for the MMT spectrum as our best spectroscopic esti¬ 
mates of the surface parameters of PG 1219+534. To repeat, 
those are = 33,600 ± 370 K, logg = 5.810 ± 0.046, 
and log A(He)/A(H) = -1.49 ± 0.08. These values are con¬ 
sistent with those derived by Koen et al. (1999). Compared to 
the Heber et al. (2000) results, our values indicate a surface 
gravity somewhat lower than proposed by them. Nonetheless, 
we point out that a much closer agreement exists with their 
“NLTE: H+He” measurement (see their Table 5). In fact, for 
the reasons given above concerning the question of mode den¬ 
sity, we may consider that measurement as “their” best solution 
as well. Einally, we stress that the derived atmospheric parame¬ 
ters are typical of most EC14026 stars and place PG 1219+534 
close to the center of the theoretical instability strip discussed 
by Charpinet et al. (2001). 

2.2. Time Series Photometry 

We observed PG 1219+534 in optical “white light” fast pho¬ 
tometric mode at the CEHT during a dedicated four night run 
scheduled on 2004, March. Table 1 provides a journal of these 
observations. As with our previous programs, the photometric 
data were gathered with Lapoune, the portable Montreal three- 
channel photometer. The data were then carefully corrected 
for various undesirable effects (e.g., atmospheric transparency 
variations, extinction, etc.) using the data reduction package 


Oscar specifically designed by one of us (P.B.) for that pur¬ 
pose. We refer the interested reader to Billeres et al. (1997) for 
details on the observations and data reduction procedures. 

The observing conditions varied from relatively poor for 
the first two nights to simply outstanding for the rest of the 
run. The quality of the data gathered during these observations 
is illustrated on a night-by-night basis in Eigure 3, where the 
left panel shows the fully reduced light curves, and the right 
panel displays the corresponding Eourier amplitude spectra in 
the relevant 0-20 mHz frequency bandpass. The light curves 
from the first and second nights (runs cfh-104, cfh-105a and 
cfh-105b), as well as from the very end of the third night (run 
cfh-106b) were rather noisy, and this is also clearly reflected 
in their associated Eourier amplitude spectra. These time series 
were obtained during bright time (contrary to our request) and 
the presence of dust and high winds as well as thin cirrus during 
the first two nights is responsible for the much lower quality of 
the data we obtained. Our past experience at the CEHT shows 
that photometric data gathered with Lapoune can be quite good 
even in presence of thin cirrus, but this works only under dark 
sky conditions. The contrast is particularly sharp when com¬ 
paring these data with the very high quality photometry (more 
typical of what we have been accustomed to at CEHT) acquired 
during most of the third night (run cfh-106a) and during the 
fourth night (run cfh-107), where conditions became photomet¬ 
ric. That last light curve in particular, a chunk of 9.36 h with no 
interruption obtained during perfectly photometric conditions, 
is among the best that have been obtained with Lapoune at 
CEHT, so far. It is shown in greater detail in Eigure 4. The vari¬ 
ations have typical peak-to-peak amplitudes of less than ~ 6% 
and are dominated by a pseudo-period of ~ 130 s. It is obvi¬ 
ous however that several other periodicities are involved in the 
variations, considering the strong beating structures that can be 
seen in the light curve and which are typical of destructive and 
constructive interferences between modes. The multiperiodic 
nature of the variations is, of course, confirmed in the Eourier 
amplitude spectra shown in the right panel of Eigure 3, where 
four dominant peaks clearly emerge over the mean noise level. 
Our thorough analysis of the light curves described in the fol¬ 
lowing section will reveal that even more periodicities, indeed, 
contribute to the luminosity variations seen in PG 1219+534. 

3. Analysis of the Light Curves 

Considering that the quality of the light curves obtained on 
PG 1219+534 during our CEHT run is particularly inhomo- 
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Figure 3. Left panel. Optical (“white light”) light curves of PG 1219+534 obtained at CFHT between 2004, March 9-12 {top to bottom slices). 
The light curves are expressed in terms of residual amplitude relative to the mean brightness of the star. Right panel. Fourier amplitude spectra 
(given as a percentage of the mean brightness of the star) in the 0-20 mHz bandpass associated with each light curve shown in the left panel. 


0 1000 2000 3000 4000 



Time in seconds (Base = BJD2453076.75000) 


Figure 4. Expanded view of the best optical light curve obtained on 
PG 1219+534 during the night of 2004, March 12. Each band covers 
an interval of 4,000 s and the amplitude is expressed in terms of the 
residual relative to the mean brightness of the star. 

geneous, we carried out detailed frequency analyses on various 
subsets of the observations. We analyzed the time series in a 
standard manner by combining Fourier analysis, least-squares 
fits to the light curves, and prewhitening techniques. The pro¬ 
cedure to extract pulsation modes from the light curves is now 
greatly eased by a new dedicated software named Berthe, de¬ 
veloped recently by one of us (P.B.). This software provides 


efficient tools to conveniently proceed along the procedure de¬ 
scribed below. 

The first step of the procedure is to compute the Fourier 
transform of the light curve. We typically use 10,000 frequency 
points per chunk of width 1 mHz in the Fourier space, which 
largely oversamples that domain. The upper slices in the panels 
shown in Figure 5 and Figure 6 illustrate the Fourier spectrum 
of the light curves of PG 1219+534 in the relevant 3-16 mHz 
bandpass. As noted previously, the light curve shows multiperi- 
odic harmonic oscillations which is reflected through the pres¬ 
ence of multiple peaks in the Fourier spectrum. We next isolate 
from the Fourier spectrum the frequency (period) of the usually 
(but not always) largest peak seen in a given frequency inter¬ 
val. That frequency (period) is used in a least-squares routine 
that provides the amplitude and the phase of the harmonic os¬ 
cillation of that given period which best-fits the complete light 
curve. This harmonic oscillation is then subtracted from the 
light curve in the time domain (prewhitening) and the Fourier 
transform of the residual light curve is computed. A second 
frequency, again often associated with the highest peak in the 
Fourier transform of the residual, is identified and the least 
squares routine is reapplied but this time with the two fre¬ 
quencies, amplitudes, and phases fitted simultaneously. This 
provides upgraded estimates for the amplitude and phase of 
the first harmonic oscillation and initial estimates for the sec¬ 
ond one. These two harmonic oscillations are then subtracted 
from the original light curve through prewhitening and another 
Fourier transform of the residual light curve is computed. The 
procedure is repeated until no useful information can be ex¬ 
tracted from the residual light curve through this method. It is 
standard to stop the iteration when the amplitude of the candi¬ 
date peaks in the Fourier domain gets below 3.7 times the local 
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rms noise level, which corresponds to the 99% statistical sig¬ 
nificance limit. We find, in practice, that applications of more 
sophisticated statistical methods to test for the reality of peaks 
in the Fourier domain (such as the false-alarm probability for¬ 
malism) generally confirm that an amplitude 3 times larger than 
the mean noise level is sufficient to account for the reality of a 
mode. 

We first analyzed the light curves obtained on 
PG 1219-1-534 separately, on a night-by-night basis, be¬ 
fore combining them according to their respective mean noise 
level, i.e., night 1 was combined with night 2, and night 3 


(without the noisy chunk at the end of the night that would 
just ruin the signal-to-noise ratio, if added) was combined with 
night 4. Finally, we examined the complete light curve built 
from all nights. For the various cases considered, illustrations 
of the Fourier transforms obtained before and after prewhiten¬ 
ing of the identified frequencies are provided in Figure 5 
and Figure 6. The periods extracted from these analyses are 
listed in Table 2 (each identified period is labeled fj, where 
j is ordered by decreasing amplitude), as well as the results 
obtained previously by Koen et al. (1999). In this table, we 
also provide information relative to the run considered. This 
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Figure 6. Same as Figure 5, but for the combined time series from night 3 (run cfh-106a only) and night 4 {left panel), and for all light curves 
assembled together {right panel). 


Table 2. Periods of the Modulations Detected in the Light Curves of PG 1219+534. 



Koen et al. 

Night 

Night 

Night 

Night 

Nights 

Nights 

Nights 


1999 

1 

2 

3 

4 

1 &2 

3&4 

All 

Resolution (pHz) 

3.64 

45.02 

35.03 

62.11 

29.68 

8.85 

8.61 

3.44 

Duty cycle (%) 

17.8 

100 

98.3 

100 

100 

44.5 

42.9 

34.4 

FT Noise (%) 

0.040 

0.040 

0.055 

0.0090 

0.0040 

0.040 

0.0055 

0.025 


Period 

Period 

Period 

Period 

Period 

Period 

Period 

Period 

Id. 

(s) 

(s) 

(s) 

(s) 

(s) 

(s) 

(s) 

(s) 

fi 




172.241 

172.164 


172.214 


/9 




159.046 

158.793 


158.789 


/4 

148.776 

148.603 

148.831 

148.796 

148.776 

148.769 

148.775 

148.777 

ft 

143.649 

143.619 

143.585 

143.639 

143.650 

143.652 

143.649 

143.647 

h 

[134.960] 



134.962 

135.147 


135.160 


h 

133.510 

133.497 

133.490 

133.519 

133.517 

133.503 

133.516 

133.518 

k 




129.239 

129.130 


129.099 


fi 

128.078 

128.115 

128.104 

128.069 

128.079 

128.075 

128.077 

128.080 

f5 




122.404 

122.397 


122.408 


[/lo] 





[82.292] 


[82.261] 


[/ll] 





[71.865] 


[71.879] 


[/ 12 ] 





[69.193] 


[69.204] 



includes the resolution achieved (in /iHz), the rms noise level 
reached in the data (in % of the mean brightness of the star), 
and the duty cycle (in % of the total coverage). Of course, 
each night considered individually gives a duty cycle of 100% 
(except for the second night, due to a short interruption in 
the time series), which then decreases due to the daily gaps 
that are introduced when data sets of successive nights are 
combined (down to 34.4% for the whole run). In the meantime, 
however, the frequency resolution improves from 29.7 /rHz 
for our longest contiguous light curve (night 4) down to 3.44 
pHz with the time baseline of the whole photometric run. We 


note that the best frequency resolution reached with our data 
set is comparable to the resolution obtained in Koen et al. 
(1999), but our coverage is better. Unfortunately, no significant 
improvement of the mean noise level in the Fourier space can 
be achieved with the whole data set compared to Koen et al. 
(1999) data. This is a direct consequence of the heterogeneous 
quality of the CFHT light curves. The rms noise level in the 
Fourier domain is approximately 5-10 times higher for the 
first two nights compared to the excellent data obtained during 
nights 3 and 4. Mixing all time series to construct the global 
light curve significantly degrades the S/N ratio in the Fourier 
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domain, as the mean noise level is essentially regulated by the 
light curve of worst quality. Obviously, in the present case, 
focussing on the best data only (i.e., those obtained during 
night 3 and night 4) would yield the best results for extracting 
pulsation modes fromPG 1219+534 photometric variations. 

The light curve from the fourth night is by far the best 
obtained during these observations and thus, not surprisingly, 
gives the Fourier transform with the lowest rms noise level 
(~ 0.004%). From that time series alone, 9 frequencies (peri¬ 
ods) could be extracted without any ambiguity. The prewhiten¬ 
ing sequence for this particular light curve is illustrated in the 
successive bands shown in the lower-right panel of Figure 5. 
The upper band displays the full amplitude spectrum in the rel¬ 
evant 3-16 mHz frequency bandpass where photometric activ¬ 
ity is concentrated. We note that from 16 mHz up to the Nyquist 
frequency (50 mHz, as we used a sampling time of 10 s), the 
structures found in the Fourier amplitude spectrum are entirely 
consistent with noise. At the low frequency end, the rms noise 
level increases due to residual atmospheric variations (which 
are difficult to remove completely) and no obvious periodici¬ 
ties could be identified (but see below). Four dominant periods, 
labeled /i to / 4 , are clearly seen at that stage, while the noise 
level is so low at the provided scale that it is essentially invis¬ 
ible. The mid-band shows the Fourier transform of the resid¬ 
ual light curve after prewhitening of those four main harmonic 
oscillations. This leads to the unambiguous detection of five 
additional periods to /g) that clearly emerge above the de¬ 
tection threshold of 3.7 times the rms noise level materialized 
by the horizontal dashed line (the dotted line indicates the limit 
of 3 times the rms noise level). The lower band displays the 
Fourier transform of the residual light curve after prewhitening 
of the nine identified periodicities. Little power is then left in 
the Fourier transform of the residual light curve, but hints of 
remaining structures exist, in particular at the high-frequency 
end of the bandpass considered. Three more structures (/lo, /i i, 
and fix) are indeed found slightly above, or just below the de¬ 
tection threshold (dashed line). These are also visible in the 
Fourier transform of the residual light curve of the combined 
nights 3 and 4 (lower-band in left-panel of Figure 6), but un¬ 
fortunately cannot be confirmed with the analysis of the other 
light curves, due to their much higher mean noise level. For this 
reason, we consider as marginal the detections of these periods, 
although we will provide below arguments which suggest that 
these structures may be real. These marginal period detections 
are indicated within brackets in Table 2. 

Globally, we find that the four dominant modes, /i to / 4 , 
are detected in all the CFHT light curves (and combinations of 
light curves) and correspond to the four periods already iden¬ 
tified in Koen et al. (1999). These authors had also suspicions 
concerning a fifth period at 134.96 s (that corresponds to /g, 
here), which we indeed confirm with no ambiguity on the ba¬ 
sis of the best light curves we obtained. Note that only the four 
dominant periods are seen when nights 1 and 2 are considered 
either separately or combined into larger data sets, such as the 
complete light curve built from all nights. This is due, again, 
to the relatively poor quality of these specific data and the dra¬ 
matic impact it has on the mean noise level in the Fourier do¬ 
main. Nonetheless, the analysis of the best available time se¬ 


ries yields a significant improvement in terms of the number of 
modes detected compared to the original data, since five addi¬ 
tional periods (A to fg) can be securely identified in the spec¬ 
trum ofPG 1219+534. 

We provide in Table 3 a list of the detailed properties that 
characterize the nine harmonic oscillations uncovered. These 
are the frequency (in /rHz), the period (in second), the am¬ 
plitude (in % of the mean brightness), and the phase (in sec¬ 
onds). The three marginal detections mentioned above are also 
included within brackets. We note that two of these can be 
interpreted as harmonics or cross-frequencies of the domi¬ 
nant peaks, i.e., /n = 2/2 and fix - fx + f 3 (and therefore 
would not constitute independent modes usable for asteroseis¬ 
mology), while the third period (/lo) may be an independent 
pulsation mode. We have, however, explicitly excluded this pe¬ 
riod and considered only the nine well secured pulsation modes 
for the following asteroseismic analysis of PG 1219+534. We 
have adopted the list of periods derived from the analysis of 
the combined light curves from night 3 and night 4, which al¬ 
lows to recover all the harmonic oscillations seen in the data 
taken during night 4 treated separately (the best achievement 
in terms of S/N ratio) while reaching a decent frequency reso¬ 
lution of 8.61 juHz. The measured frequencies (each associated 
with a frequency peak in the Fourier domain) are probably ac¬ 
curate to 1/10 of the formal resolution (itself corresponding to 
the width at half-maximum of such a peak), i.e., to within 0.86 
juHz (see Bretthorst 1988). The amplitudes and phases derived 
through the least-squares technique have formal errors as given 
in Table 3. We note at this point, that the gain in frequency reso¬ 
lution obtained while combining all light curves for the Fourier 
analysis does not reveal the presence of close frequencies near 
the dominant peaks (due to rotational splitting, for instance). 
This observation holds for the five newly discovered frequen¬ 
cies as well, although the frequency resolution is lower in these 
cases. Such a finding indicates that PG 1219+534 is likely a 
slow rotator with a rotation period longer than ~ 3.4 days, i.e., 
the longest time baseline achieved with the photometry cur¬ 
rently available. Finally, we note a significant change in the 
amplitude of three of the four dominant periods that have been 
detected in both our contemporary data and the older Koen 
et al. (1999) light curves. Except for the period fi which is 
found slightly weaker (by ~ 6%) in our time series, the periods 
fx, f}, and /4 all appear much stronger than observed in Koen 
et al. (1999), by factors up to ~ 2-4. Such amplitude vari¬ 
ations are not uncommon among rapidly pulsating sdB stars. 
In PG 1219+534, these may be due to beating between very 
close - thus, unresolved with current data - frequencies (multi- 
plet components caused by slow rotation, for instance), and/or 
to significant changes in the intrinsic amplitudes of the modes 
themselves. However, additional observations with longer time 
baselines will be needed to further discuss the stability of pul¬ 
sation mode amplitudes in this star. 

Interestingly, the odd conditions under which the CFHT 
data were obtained provides us with a perfect illustration of 
the fact that mixing low and high S/N ratio time series can es¬ 
sentially ruin all the benefit of having the high-sensitivity data. 
The overall mean noise level in the Fourier space is clearly 
mostly regulated by the noise level of the worst light curve in- 
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Table 3. Properties of the Harmonic Oscillations Identified in PG 1219+534. 


Id. 

Frequency 

ffiHz) 

Period 

(s) 

Amplitude 

(%) 

Phase 

(s) 

Comments 


fi 

5806.71 

172.214 

0.0456 ± 0.0030 

68.72+ 1.80 



f9 

6297.65 

158.789 

0.0334 ± 0.0030 

23.14 + 2.27 



/4 

6721.55 

148.775 

0.4053 ± 0.0030 

42.86 + 0.18 

detected with an amplitude of 0.22% by Koen et al. 

(1999) 

fl 

6961.39 

143.649 

0.7228 ± 0.0030 

102.50 + 0.09 

detected with an amplitude of 0.18% by Koen et al. 

(1999) 

h 

7398.62 

135.160 

0.0429 ± 0.0030 

98.73 + 1.51 



h 

7489.73 

133.516 

0.6874 ± 0.0030 

43.65 + 0.09 

detected with an amplitude of 0.23% by Koen et al. 

(1999) 

U 

7745.98 

129.099 

0.0544 ± 0.0030 

104.69+ 1.13 



fl 

7807.79 

128.077 

0.7972 ± 0.0030 

44.07 + 0.08 

detected with an amplitude of 0.85% by Koen et al. 

(1999) 

h 

8169.43 

122.408 

0.0963 ± 0.0030 

95.70 + 0.61 



[/lo] 

[12156.40] 

[82.261] 

[0.0190 ± 0.0030] 

[59.38 + 2.07] 

slightly above detection threshold; not used for seismology 

Uni 

[13912.31] 

[71.879] 

[0.0143 ± 0.0030] 

[21.06 + 2.40] 

2 / 2 ; not used for seismology 


Uni 

[14450.08] 

[69.204] 

[0.0188 + 0.0030] 

[48.10+ 1.76] 

fl + h \ not used for seismology 



eluded in the analysis. It is instructive that, for this particular 
star PG 1219+534 which, at the outset, turns out to have a rel¬ 
atively simple and well resolved pulsation spectrum, just one 
night of high sensitivity CFHT data appears to be sufficient to 
make this object amenable to detailed asteroseismic analysis. 
More generally, for the purpose of asteroseismology of stars 
with simple enough oscillation spectra (such as EC 14026 stars 
often have, with some exceptions), high sensitivity is often to 
be sought preferrentially to frequency resolution and/or cov¬ 
erage. This, indeed, has some implications for the treatment 
of multisite campaign data where light curves of various qual¬ 
ity are generally mixed. In those cases, mixing low S/N data 
with high S/N data, although improving resolution and reduc¬ 
ing aliasing problems, can significantly degrade the S/N ratio 
reached with the best available data. Accordingly, weighting 
techniques may contribute, in some circumstances, to reduce 
the deterioration of the overall S/N ratio, although generally at 
the expense of a degraded window function (see, for instance, 
the discussion of Handler 2003). However, with very heteroge¬ 
neous data sets, these techniques would still lead to a significant 
increase of the mean noise level in the Fourier domain com¬ 
pared to the S/N ratio achieved with the best light curves. This 
has to be considered carefully when the goal is to extract low 
amplitude variations to increase the number of modes for aster¬ 
oseismology. Of course, the best is to have both high sensitivity 
and coverage with data as homogeneous as possible. Along this 
line, we hope for future bi- (or tri-) site fast photometric obser¬ 
vations based on 4m-class telescopes, as it would constitute, in 
addition to other developing techniques (multicolor photome¬ 
try and high time resolution spectroscopy), a significant step 
forward in the field of asteroseismology of EC14026 stars. 

To end this Section, let us briefly mention that, while per¬ 
forming the Fourier analysis of the best light curves obtained 
on PG 1219+534, we have discovered hints of low frequency 
modulations in the amplitude spectrum. Having a closer look 
at the highest S/N ratio Fourier amplitude spectra shown in the 
right panel of Figure 3 (especially the lowest band correspond¬ 
ing to the data of our best night), one can notice low amplitude 
structures below 1 mHz that still significantly emerge above 
the local mean noise level. If such periodic modulations are 


real and due to pulsations, those peaks would correspond to 
modulations with periods much longer than the typical peri¬ 
ods observed in EC 14026 stars, thus involving relatively high- 
order gravity modes. To some extent, it would be similar to 
the variations seen in the long period sdB pulsators, although 
PG 1219+534 does not have, in principle, the expected physi¬ 
cal parameters to belong to that class of pulsating stars which, 
typically, have much lower effective temperatures and surface 
gravities. A possible analogy would be with the recent excit¬ 
ing discovery of a long-period photometric variation in the 
EC 14026 star HS 0702+6043 (Schuh et al. 2005) which sug¬ 
gests that this star may pulsate in both low-order p-modes and 
high-order g-modes. It is, however, premature to speculate fur¬ 
ther on the nature of these (potential) low frequency luminos¬ 
ity modulations in PG 1219+534. Additional photometric ob¬ 
servations at high sensitivity with a longer time baseline than 
is currently available will be necessary to address this partic¬ 
ular issue and possibly confirm (or rule out) the presence of 
g-modes in PG 1219+534. Specific efforts along this line are 
underway. 

4. Asteroseismic interpretation of the 
Observations 

4.1. The Double-Optimization Scheme 

The asteroseismological interpretation of the period spectrum 
of pulsating stars has often been impaired by the lack of identi¬ 
fication of the modes being observed. The analysis of EC 14026 
stars is confronted to similar challenges as white light fast pho¬ 
tometry generally gives little clue - and quite often no clue at 
all - about the degree i and/or radial order k of the modes de¬ 
tected. Indeed, for PG 1219+534, we have no a priori knowl¬ 
edge of the geometry of the nine observed modes (not even 
hints on the i value, since no splitting due to rotation has been 
detected at the current frequency resolution). Nonetheless, we 
have pursued major efforts in the last few years to develop a 
new method that allows us to bypass this difficulty. 

Our approach is based on the well known forward method 
which consists of comparing directly periods computed from 
stellar models to periods observed in the star under considera- 
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tion with the goal of reproducing as accurately as possible its 
oscillation spectrum. Applied to EC 14026 stars, it consists of a 
double-optimization procedure that takes place simultaneously 
at the period matching level and in the model parameter space. 
Due to the lack of mode identification from the available data, 
assigning No observed periods to Nt theoretical periods from a 
given model {Nt > No, in general) usually leaves a myriad of 
possible combinations. The first optimization therefore consists 
of finding the combination (or mode identification) that leads to 
the best possible simultaneous match of all the observed peri¬ 
ods for that given model. We note that this best match does not 
necessarily have to be a good match at this point. This will be 
the role of the second optimization procedure discussed below. 
More formally, the quality (or merit) of a period fit is evaluated 
as the sum of the squared differences between the observed 
periods and their assigned theoretical periods (i.e., through a 

least-squares formalism). The chosen mode identification for a 

_2 

given model is then the combination that provides a x value 
corresponding to the minimum of that sum among all the pos¬ 
sible associations, that is 


N„ 


X =mm<x 


= Z 


d(0 

obs 


d(0 \ 


, V combinations|. (1) 


Here, is the observed period and represents its 
associated theoretical period, cr, is an optional weight that 
can be attributed either individually on periods or globally 
(see Subsection 4.3). The second optimization is carried out 
in the model parameter space where the quantity x derived 
from the first optimization is now seen as a function of the 
N model parameters (considering here the general case), i.e., 
X = X {oi,a 2 , ■■■,0^), where a, represents one of the N pa¬ 
rameters. The procedure then consists of localizing the opti¬ 
mal model (or models) that minimizes this “merit” function in 
the A-dimensional space, leading objectively to the best period 
match of the observations. 

We have developed various sets of numerical tools that al¬ 
lows us to apply this double-optimization scheme to isolate best 
period fitting models for asteroseismology. Our original efforts 
along those lines were described in Brassard et al. (2001). Since 
then, many important improvements have been made and two 
independent packages (fortunately leading to the same results!) 
are now operational, one in Montreal and the other in Toulouse. 
Detailed descriptions of the numerical tools developed in these 
packages will be reported elsewhere. Briefly, the Toulouse 
package extensively used for this analysis of PG 1219+534 
includes a genetic algorithm based period matching code that 
performs the first optimization part, a parallel grid computing 
code useful to explore the parameter space and visualize the 
complex shape of the x (oi, 02 ,ov) function, and a paral¬ 
lel multimodal genetic algorithm based optimization code de¬ 
signed to explore efficiently the vast A-dimensional parameter 
space and seek for the optimal solutions. By locating simul¬ 
taneously the global and local minima of the x^ function, this 
code is capable of identifying several optimal models if more 
than one solution exists. 

It is important to realize that this approach is in fact built on 
the very simple requirement that models pretending to provide 


a good asteroseismological fit of the spectrum of a pulsating 
star must match all the observed periods simultaneously, i.e., 
the procedure is a global optimization. Moreover, we stress that 
the solutions identified from this approach do not rely on any 
previous mode identification. The mode identification appears 
instead naturally as the solution of the global fitting procedure, 
i.e., the mode identification obtained is the one that provides 
the best simultaneous match of all the detected periods. 


4.2. Equilibrium Models and Computations of their 
Pulsation Properties 


Applications of our double-optimization scheme to EC14026 
stars requires computing appropriate theoretical period spec¬ 
tra that need to be compared with the observed periods. Three 
codes are involved in this process, starting with an equilibrium 
model building code first described in Brassard & Eontaine 
(1994) and adapted to produce the so-called “second gener¬ 
ation” models suitable for pulsating sdB stars (see Charpinet 
et al. 1997, 2001). These models are improved static envelope 
structures extending as deep as log^ = log[l - M(r)/M,] ^ 
-0.05. Eor the purpose of pulsation calculations, the central 
missing nucleus that contains ~ 10% of the stellar mass is 
treated as a “hard ball”. This approach, in comparisons to 
full evolutionary stellar models, gives excellent pulsation re¬ 
sults, especially for the relatively shallow p-modes that are, 
in hot B subdwarf stars, insensitive to the detailed structure 
of the central parts of the core as was discussed at length in 
Charpinet (1999) and Charpinet et al. (2000, 2002). Moreover, 
these “second generation” static models are superior to all 
evolutionary structures currently available in that they incor¬ 
porate the nonuniform iron abundance profiles predicted by 
the theory of microscopic diffusion assuming an equilibrium 
between gravitational settling and radiative levitation. We re¬ 
call that diffusion leads to the constitution of a reservoir of 
iron in the H-rich envelope of subdwarf B stars with signif¬ 
icant iron enrichments produced locally which are responsi¬ 
ble, through the /e-mechanism, for the destabilization of the 
low-order, low-degree p-modes observed in EC 14026 pulsators 
(Charpinet et al. 1997). This additional ingredient in the con¬ 
stitutive physics of the models is therefore essential to repro¬ 
duce the excitation of the pulsation modes. In addition, micro¬ 
scopic diffusion modifies sufficiently the stellar structure to in¬ 
duce significant changes to the pulsation periods themselves, 
thus impacting on the asteroseismic analysis. Consequently, 
radiative levitation is a crucial ingredient that must be taken 
into account in the modeling of pulsating sdB stars for the pur¬ 
pose of accurate asteroseismology. Eour fundamental parame¬ 
ters are required to specify the internal structure of hot B sub¬ 
dwarf stars with the second-generation models: the effective 
temperature Teff, the surface gravity logg (traditionally given 
in terms of its logarithm), the total mass of the star M,, and 
the logarithmic fractional mass of the hydrogen-rich envelope 
log^(H) s log[M(H)/M,]. The latter parameter is intimately 
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related to the more familiar parameter Menv, which corresponds 
to the total mass of the H-rich envelope'. 

In a second step, we use two efficient codes based on fi¬ 
nite element techniques to compute the pulsation properties of 
the model. The first one is an updated version of the adiabatic 
code described in detail in Brassard et al. (1992). It is used 
as an intermediate (and necessary) step to obtain estimates for 
the pulsation mode properties (periods and eigenfunctions) that 
are then used as first guesses in the solution of the full, nona- 
diabatic set of nonradial oscillation equations. The second one 
is an improved version of the nonadiabatic code that has been 
described briefly in Fontaine et al. (1994) and which solves 
the nonadiabatic eigenvalue problem. It provides the necessary 
quantities to compare with the observations. For the purpose of 
asteroseismology, these are mainly the periods and the stability 
coefficients. To illustrate the type of quantities that are derived 
from the whole process, we provide, in Table 4, a typical out¬ 
put of the nonadiabatic pulsation code that gives the pulsation 
properties of a given model. While it refers specifically to the 
optimal model (see the discussion in the next subsection), we 
only wish to emphasize the illustrative aspect of the theoret¬ 
ical results at this point. For each mode found in the chosen 
period interval. Table 4 gives the degree (, the radial order k, 
the period Pth (= where ctr is the real part of the com¬ 

plex eigenfrequency), the stability coefficient ct/ (the imaginary 
part of the eigenfrequency), the logarithm of the so-called ki¬ 
netic energy of the mode E, and the dimensionless first-order 
solid rotation coefficient Cm- As is standard, our equilibrium 
stellar models are perfectly spherical, and each mode is 21 + 1 
fold-degenerate in eigenfrequency. 

For asteroseismic studies, the most important quantities de¬ 
rived from theoretical calculations are of course the periods of 
the modes. These are sensitive to the global structural param¬ 
eters of a model which we seek to infer for a real star through 
a comparison with a set of observed periods. Because nona¬ 
diabatic effects on the periods are small (but included in our 
calculations), asteroseismology could very well be carried out 
at the level of the adiabatic approximation only. However, the 
stability coefficient, a purely nonadiabatic quantity, also pro¬ 
vides useful information on the mode driving mechanism. A 
positive value of crj indicates that a mode is damped (or “sta¬ 
ble”), and therefore that it should not be observable. A negative 
value of CT/ indicates that the mode is excited (or “unstable”) in 
the model, and consequently that it may reach observable am¬ 
plitudes. As illustrated in Table 4, modes are excited within a 
band of periods associated with low-order acoustic modes. This 
is typical of all our models of rapid sdB pulsators. In terms of 
asteroseismology, it is of course within that bandpass that one 
would expect to find the observed periods of a real EC14026 
star and this important stability information justifies the use of 


* We stress that the parameter commonly used in Extreme 
Horizontal Branch stellar evolution includes the mass of the hydrogen 
contained in the thin He/H transition zone, while the parameter M(H) 
in our models does not. They are related through log[Menv/M,] = 
log[M(H)/M,] + C, where C is a small positive term, slightly depen¬ 
dent on the model parameters, determined by the mass of hydrogen 
that is present inside the He/H transition zone itself. 


full nonadiabatic calculations in our case. For its part, the ki¬ 
netic energy is a secondary quantity in the present context. It 
provides a measure of the required energy to excite a mode 
of a given amplitude at the surface of a star. Since this is a 
normalized quantity, only its relative amplitude from mode to 
mode is of interest. The kinetic energy is a useful diagnostic 
tool for determining where pulsation modes are formed. Larger 
values of E imply that such modes probe deeper into the star, 
in higher density regions. The kinetic energy bears the signa¬ 
ture of mode trapping/confinement, which are caused by par¬ 
tial reflexions of propagating waves onto thin chemical transi¬ 
tion zones (between the pure He-core and the H-rich envelope, 
mainly). Finally, the rotation coefficient is useful for interpret¬ 
ing the fine structure in the Fourier domain in terms of slow 
rotation that lifts the 2£ + I fold-degeneracy of the mode peri¬ 
ods. At the frequency resolution reached in our observations of 
PG 1219+534, no such structures have been detected, however. 

4.3. The Quest for the Optimal Model of PG 1219+534 

The search for the optimal model solution (or solutions) that 
can best reproduce the nine identified periods of PG 1219+534 
was initiated with our multimodal genetic algorithm (GA) opti¬ 
mization code. As mentioned briefly in the previous subsection, 
the search is carried out in the four-dimensional space defined 
by the main parameters T^g, logg, M,, and log^(H) that spec¬ 
ify the structure of a sdB star with our second-generation mod¬ 
els. The GA code, designed for efficient explorations of a wide 
parameter space domain, was first used to localize potential re¬ 
gions were solutions may exist. For that purpose, initial bound¬ 
aries of the search domain were defined as follows : 37,000 K 
> Teff > 28,000 K, 5.65 < logg < 5.95 for the surface param¬ 
eters, and -5.0 < log^(H) < -2.0, 0.30 < M,/Mo < 0.53 for 
the structural parameters. 

These limits were loosely set according to current spectro¬ 
scopic estimates for the atmospheric parameters Tefr and log g 
of PG 1219+534. We emphasize the fact that any eventual de¬ 
termination of the atmospheric parameters of PG 1219+534 
through asteroseismic means has to be consistent with our 
spectroscopic estimates. Since the latter are based on standard, 
well-understood techniques (see, however, the discussions of 
Wesemael et al. 1997 and Heber & Edelmann 2004 on possible 
systematic errors introduced by these methods when applied to 
sdB stars), the validity of our pulsating models would have to 
be questioned if the results show otherwise. Constraints on the 
two other parameters M, and log ^(H) rely on stellar evolution 
theory. The evolutionary calculations of Dorman et al. (1993) 
indicate that hot B subdwarfs are core helium burning stars on 
the Extreme Horizontal Branch that evolve as “AGB-Manque" 
objects. According to these calculations, their possible masses 
are found in a narrow range 0.40 - 0.43 < M,/Mq < 0.53, 
with a somewhat uncertain lower limit related to the mini¬ 
mum mass required to ignite helium and a more sharply de¬ 
fined upper limit above which the models evolve to the AGB. 
More recently, however, Han et al. (2003) suggested, in their 
investigation of various binary evolution scenarios for the for¬ 
mation of sdB stars, a somewhat larger distribution of masses 
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among these stars. Although these authors derived mass distri¬ 
butions which are strongly peaked around the canonical value 
of ~ 0.47 Mq, they suggest that the mass of an sdB star could 
be as low as ~ 0.30 Mq if formed through the Roche lobe over¬ 
flow channel or as high as ~ 0.70 Mq if formed through the 
merger channel. To explore these possibilities with asteroseis¬ 
mology, we adopted a generous lower mass limit of 0.30 Mg 
but kept 0.53 Mg as our upper boundary (results from test cal¬ 
culations for models with higher masses will be briefly dis¬ 
cussed, though). Finally, the range of possible values for the 
last parameter log ^(H), which is related to the mass of the FI- 
rich envelope, was chosen according to the work of Dorman 
et al. (1993) in order to fully map the region of the logg - Teff 
plane where sdB stars are found. 

For the pulsation calculations step, we considered all the 
modes (be it of p-, /-, or g-type) of degree f = 0 up to = 3 
with periods in the range 70 - 500 s, i.e., covering amply the 
range of periods observed in PG 1219+534. The usual argu¬ 
ment to limit the value of { above which modes are no longer 
observable is geometric cancellation on the visible disk of the 
star. A quantitative expression of the visibility argument used 
to restrict the search to low values of ( is obtained by evalu¬ 
ating numerically the geometric factor given, for example, by 
equation (B7) of Brassard et al. (1995). Using an Eddington 
limb darkening law in that formulation gives visibility factors 
with values of 1.000, 0.7083, 0.3250, 0.0625, 0.0208, 0.0078, 
0.0078, 0.0023, 0.0039, 0.0009, 0.0023, 0.0005, and 0.0015, 
respectively, for f = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, and 
12. This cancellation effect naturally favors a priori values of 
{ - 0, 1, and 2 in pulsating stars in general. We note, however, 
that the mode densities and period distributions seen in sev¬ 
eral well-observed EC14026 stars such as KPD 1930+2752, 
PG 1605+072, PG 1047+003, and PG 0014+067 force us to 
consider modes with higher £ values {{ - 3, or even {-Am 
some cases). Otherwise, there would be not enough theoretical 
modes available in the observed period window to account for 
the observations. To justify further the (small) heresy that con¬ 
sists of considering modes of degree up to ^ = 3 or 4 (although, 
again, observational evidence already suggests that such modes 
are effectively seen), we stress that our growing observational 
experience with EC 14026 pulsators indicates clearly that in¬ 
trinsic mode amplitudes in these stars are not the same from 
mode to mode. In addition, these amplitudes may vary signif¬ 
icantly on a monthly or yearly basis, some modes being visi¬ 
ble one season but disappearing the next one. Therefore, one 
should be cautious with the simple geometric cancellation ar¬ 
gument which is based on the implicit assumption that intrin¬ 
sic amplitudes of all modes are the same. We suggest instead, 
at least for EC 14026 stars, that the limit in ( at which modes 
cease to be observable due to cancellation effects cannot be a 
strict threshold, since it also depends on the intrinsic amplitude 
of the modes considered individually. Eor the purpose of aster¬ 
oseismology, we choose to consider modes from f = 0 up to the 
minimum value of the degree ( needed to account for the mode 
density in the observed period range. Eor the PG 1219+534 
periods identified in Table 3, modes up to f = 3 are needed. 
Of course, we cannot rule out completely the presence of the 
odd { - A modes (or perhaps even a mode of higher degree) 


with an unusually high intrinsic amplitude in the light curve of 
PG 1219+534. Our approach explicitly excludes this possibil¬ 
ity, however, and we seek to interpret the period distribution in 
that star solely in terms of f = 0 - 3 modes. 

Within the search domain specified, the GA code identified 
two families of solutions that best match the observed periods. 
A third family of potential solutions was also spotted near the 
low-gravity edge of the search domain, but was immediately 
rejected as being highly inconsistent with the available spec¬ 
troscopic measurements (see discussion below). These model 
solutions are essentially equivalent in terms of quality of fit, 
indicating that best-fitting the observed periods only leads to 
a degenerate situation (there is another source of solution de¬ 
generacy which will be discussed below). The optimal model 
corresponding to the first family of solutions was found at 
Teff = 33,640 K, logg = 5.8071, log^(H) = -4.3448, and 
M, = 0.4570 Mq with a value of 0.0784. The model as¬ 
sociated with the second family of solutions was localized at 
Teff = 33,715 K, logg = 5.7178, log^(H) = -3.8142, and 
M, = 0.4680 Mq with a)^^ value of 0.0878. Note that in the 
evaluation ofx , we have chosen a global weight factor cr, = cr, 
where cr is the inverse of the theoretical mode density, i.e., the 
ratio of the width of the period window (here 500 -70 - 430 s) 
to the number of theoretical modes in that window for a given 
model. This choice was made in order to ease the automatic 
search of the best-fit models in the parameter space by offset¬ 
ting, at least in part, the built-in bias in favor of models with 
a higher theoretical mode density (generally, the models hav¬ 
ing a lower surface gravity). Of course, this choice does not 
influence the positions of the minimum values of in the pa¬ 
rameter space domain. Maps displayed in Eigure 7 and Eigure 
8 illustrate the complex shape of the;^^ -function (shown as iso¬ 
contours of constant value of log^^) in the vicinity of the two 
potential solutions, whose exact locations according to the GA 
code are indicated by black crosses. These figures respectively 
show slices of this function along the logg - Teg plane (at fixed 
parameters M* and log< 7 (H) set to their optimal values) and 
along the M, - log ^(H) plane (at fixed parameters and log g 
set to their optimal values). Best fitting models corresponding 
to low values of;^f appear as dark blue regions, while red areas 
indicate regions of the model parameter space with high values 
ofx , i.e., where theoretical periods computed from models do 
not fit well the observed periods. Considering the logarithmic 
scale used to represent the merit function on these plots, we 
point out that the blue regions correspond to well-defined min¬ 
ima. 

If one family of solutions among the two identified cannot 

_2 

be preferred over the other one simply on the basis of the x 
value alone, additional constraints provided by spectroscopic 
measurements of the atmospheric parameters of PG 1219+534 
come in handy for deciding which of these two best-fitting re¬ 
gions corresponds to the “correct” one. The values for logg 
and Teff given by Koen et al. (1999) and their associated 
errors (shown as the dashed-line rectangle in Eigure 7) do 
not permit us to clearly discriminate between the two solu¬ 
tions. However, the improved measurements obtained from 
our medium-resolution, high signal-to-noise MMT spectrum 
(shown as the solid-line rectangle in Eigure 7) allows us to un- 
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Figure 7. Left panel: Slice of the -function (in logarithmic units) along the log g - T^s plane at fixed parameters M, and log (?(H) set to 
their optimal values found for the first best-fit model solution (M, = 0.4570 Mq and log^(H) = -4.3448). The solid-line rectangle shows 
our spectroscopic estimate with its uncertainties for the atmospheric parameters of PG 1219+534, while the dashed-line and dashed-dotted- 
line rectangles represent the Koen et al. (1999) and Heber et al. (2000; from NLTE H+He model atmospheres) spectroscopic measuremen ts, 
respectively. The dotted-line rectangle shows the estimated error box associated with our asteroseismic determination of the parameters log g 
and Teff (see text). The black cross indicates the position of the best-fit solution. Right panel: Same as left panel but for the second best-fit 
model solution (the fixed parameters are set to the optimal values M, = 0.4680 Mq and logg(H) = -3.8142). 


ambiguously decide that the solution at higher surface grav¬ 
ity (i.e., the first family of solutions) is the preferred one. The 
spectroscopic values adopted for PG 1219+534 by Heber et al. 
(2000) suggests an even higher surface gravity (at logg = 
5.95 + 0.10), thus further dismissing the second family of solu¬ 
tions as a viable choice. We find, however, that no model with 
such a high surface gravity can reproduce the period spectrum 
observed in PG 1219+534. A closer look at the period spectra 
computed for high-gravity models indicates that this is sim¬ 
ply due to the fact that some of the observed periods would 
have to be associated with low-order gravity waves. However, 
the density of g-modes is too low to be compatible with the 
mode density observed in that star, and matching simultane¬ 
ously several periods with g-modes is impossible. Nonetheless, 
we again stress that the atmospheric values derived by Heber 
et al. (2000) from NLTE H+He model atmospheres (shown as 
the dashed-dotted line rectangle^ in Figure 7), which are pre¬ 
sumably the most realistic in terms of the physics involved, 
are in much closer agreement with our own spectroscopic es¬ 
timates and do not conflict with asteroseismology. In light of 
the various spectroscopic measurements, we therefore adopt 
the first identified family of models, located near logg ~ 5.81, 
as the best candidate solution that would reproduce simulta¬ 
neously the nine periods observed in PG 1219+534. We point 


^ Note that Heber et al. (2000) do not provide error estimates specif¬ 
ically associated with the parameters of PG 1219+534 derived from 
their NLTE H+He model atmospheres. We have therefore adopted 
ATeff = 500 K and A logg = 0.05 as realistic representative values. 


out at this stage that beyond the best-fit model degeneracy dis¬ 
cussed above, a second more subtle type of degeneracy of the 
asteroseismic solution among the chosen family of models does 
exist, however. This degeneracy occurs when a change in one 
of the model parameters is almost exactly compensated by a 
change in another model parameter such that the computed pe¬ 
riods remain unmodified. This phenomenon occurs in our anal¬ 
ysis of PG 1219+534 at two levels. 

First, there is a small correlation between the log ^(H) and 
Teff parameters. Starting from the derived best-fitting model, 
changing the value of log g'(H) while keeping M* constant and 
set to its optimal value generates, in the logg - Teff plane, a 
shift of the position of the local minimum which mainly fol¬ 
lows the effective temperature axis (a tiny shift also occurs 
in logg). This trend is illustrated with the map provided in 
Figure 9, which shows the “projection” of the logg'(H) axis 
onto the surface gravity-effective temperature plane. More pre¬ 
cisely, the value associated with each grid point shown on this 

map is the minimum value found among all the values of the 

_2 

X -function computed at fixed logg (set to the values as¬ 
sociated with the considered grid point), and M* (set to its op¬ 
timal value, that is M, = 0.4570 Mq) but with the parame¬ 
ter logg'(H) varying within the limits of the specified search 
domain, i.e., between logq'(H) = -5.0 and -2.0. The various 
families of solutions identified previously now appear simul- 
taneously as (blue) parallel valleys of low x in the logg-Teff 
plane. These correspond, however, to a different range of val¬ 
ues for the parameter log^(H). The small labelled axis posi- 
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Figure 8. Left panel: Slice of the ^^-function (in logarithmic units) along the M, - log qlLl) plane at fixed parameters Feff and log g set to their 
optimal values found for the first best-fit model solution (Peff = 33,640 K and logg = 5.8071). The dotted-line rectangle shows the estimated 
error box associated with our asteroseismic determination of the parameters log^(H) and M, (see text). The black cross indicates the position 
of the best-fit solution considered. Right panel: Same as left panel but for the second best-fit model solution (the fixed parameters are set to the 
optimal values T^fs - 33,715 K and logg = 5.7178). 


tioned along the valley associated with the preferred solution 
indicates the position of the local minimum of x as a func¬ 
tion of logg'(H). There is a clear, monotonic trend showing that 
this minimum shifts from higher effective temperatures (e.g., 
Teff ~ 35,800 K for log^(H) ~ -4.55) to lower T^g (~ 32,000 
K for logq'(H)-4.175) as the envelope mass of the star in¬ 

creases (i.e., logg'(H) increases). However, at the same time, 
we note a degradation of the quality of the period fit (i.e., an 
increase of the absolute value of x associated with the mini¬ 
mum) for values of log g'(H) that departs from the optimal value 
uncovered. This indicates that the degeneracy associated with 
the parameter log q'(H) has boundaries, i.e, the optimal solution 
still occupies the center of a well defined region of the parame¬ 
ter space. Interestingly, we note at this point that best-fit model 
candidates for the observed periods of PG 1219+534 do all cor¬ 
respond to models with thin H-rich envelopes. In our search, 
no acceptable asteroseismological fit has been found for mod¬ 
els with thick envelopes, i.e., with high values of the parameter 
log^(H). 

A similar map was constructed to visualize the “projec¬ 
tion” of the M,-axis onto the logg - T^g plane. It is shown 
in Figure 10. The parameter logq'(H) was kept constant (set to 
its optimal value of logg'(H) = -4.3448) and the total mass 
was varied within the limits we have defined for the search do¬ 
main, i.e., between M, - 0.53 M© and M* = 0.30 M©. The 
map indicates clearly that a correlation also exists between the 
parameter M* and the parameters T^g and, to a lesser extent, 
logg. A change in M, generates a shift in both T^g and logg 
of the position of the x^ minimum. However, contrary to the 
case of the log^(H) parameter discussed previously, we ob¬ 


serve no degradation of the quality of the period fit over the 
entire range considered for the M, parameter. Clearly in the 
case of PG 1219+534, the impact of changing the value of M, 
on the theoretical modes fitted to the observed periods can be 
completely compensated by an adequate change of both Teff 
and logg. This leads to a line-degeneracy in the ^^-function 
along which models provide comparable best-fit solutions of 
the periods. This is well illustrated in Figure 10 by the pres- 
ence of a long and flat valley of minimum value of x with 
no high-Teff and low-Teff boundaries. We also mention that a 
similar line-degeneracy associated with the second family of 
solutions, near logg ~ 5.71, does exist, although it is not ap¬ 
parent in Figure 10 as it requires a larger optimal value for the 
log^(H) parameter. The labelled axis plotted along the valley 
of minimum x^ indicates the position of the best-fit solution 
in the logg - Tefr plane for the given value of the total mass. 
The correlation between M, and {Teg, logg) is monotonic, as 
the effective temperature and surface gravity of the solution de¬ 
creases when the total mass decreases. For instance, a value of 
M, - 0.5075 M© places the minimum of^ near Teg ~ 36,600 
K and logg ~ 5.822, while a value of M, = 0.3650 M© places 
it near Teg ~ 28,400 K and logg ~ 5.796. Note that the cor¬ 
relation with the effective temperature is much stronger than 
the correlation with the surface gravity parameter. Of course, 
this trend extends beyond the limits provided in Figure 10, and 
models with lower (higher) masses lead to even cooler (hotter) 
solutions. Again, one cannot discriminate between the models 
on the basis of their ^^-value alone, and spectroscopy must be 
used to lift the degeneracy. The effective temperature derived 
from the analysis of our medium-resolution, high signal-to- 
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Figure 9. Slice of the “projected” -function (in logarithmic units) along the logg - plane at fixed parameter M, set to its optimal value 
found for the chosen best-fit model solution (M, = 0.4570 Mg). The projected, log ^(H) parameter was varied between -2.0 and -5.0 (by steps 
of 0.025). The labelled axis positioned along the valley of minimum indicates the exact location of the local minimum ofx^ for the value of 
log ^(H) discussed in the text. The solid-line rectangle shows our spectroscopic estimate with its uncertainties for the atmospheric parameters 
of PG 1219+534, while the dashed-line and dashed-dotted-line rectangles represent the Koen et al. (1999) and Heber et al. (2000; from NLTE 
H+He model atmospheres) spectroscopic measurements, respectively. 


noise MMT spectrum (the solid-line rectangle in Figure 10) al¬ 
lows us to select the appropriate section along the line of degen¬ 
eracy which corresponds to the “correct” solution. Of course, 
the necessity to rely on the spectroscopic measurement of Teff 
to uniquely derive the total mass of PG 1219+534 from as- 
teroseismological means indicates, unfortunately, that M, can¬ 
not be measured independently of T^ff for that star. Since, the 
optimal model first isolated by our GA code with parameters 
Teff = 33,640 K, logg = 5.8071, log^(H) = -4.3448, and 


M, = 0.4570 Mq is consistent with the spectroscopic determi¬ 
nation of the effective temperature of PG 1219+534, we there¬ 
fore adopt it as our preferred solution that best matches the 
period spectrum observed in this star. 

Of interest in the context of binary evolution scenar¬ 
ios to form hot B subdwarf stars, we find that the low- 
mass (high-mass) models consistent with the observed pe¬ 
riods of PG 1219+534 would require effective temperatures 
which are significantly cooler (hotter) than currently mea- 
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Figure 10. Slice of the “projected” jp^-function (in logarithmic units) along the log g - T^e plane at fixed parameter log ^(H) set to its optimal 
value found for the chosen best-fit model solution (logfjrfH) = -4.3448). The projected, M, parameter was varied between 0.30 Mq and 0.53 
Mg (by steps of 0.0025 Mg). The labelled axis positioned along the valley of minimum indicates the exact location of the local minimum of 
for the value of M. discussed in the text. The solid-line rectangle shows our spectroscopic estimate with its uncertainties for the atmospheric 
parameters of PG 1219+534, while the dashed-line and dashed-dotted-line rectangles represent the Koen et al. (1999) and Heber et al. (2000; 
from NLTE H+He model atmospheres) spectroscopic measurements, respectively. 


sured by spectroscopic means. Therefore, those solutions can 
clearly be rejected on the basis of consistency between spec¬ 
troscopy and asteroseismology. Remarkably, the mass derived 
for PG 1219+534, which is compatible with the spectroscopi¬ 
cally measured effective temperature, is found to be consistent 
with the canonical value for Extreme Horizontal Branch ob¬ 
jects. 


4.4. Period Fit and Mode Identification 

The optimal model finally isolated in the previous subsection 
provides the best simultaneous period match of the nine periods 
detected in PG 1219+534 and leads to the identification of the 
pulsation modes involved in the luminosity variations observed 
in that star. We recall that our global optimization method ap¬ 
plied to the asteroseismic analysis of PG 1219+534 does not 
rely on any previous assumption concerning the values of the 




















































18 


Charpinet et al.: Structural Parameters of the Hot Pulsating B Subdwarf PG 1219+534 from Asteroseismology 


Table 4. Pulsation Properties of the Best-Fit Model Solution and Mode Identification. 


e 

k 

Fobs 

(s) 

P,h 

(s) 

CTl 

(rad/s) 

Stability 

logs 

(erg) 

Cki 

APjP 

(%) 

Comments 

0 

6 


76.154 

+5.530 X 10-^ 

stable 

39.984 




0 

5 


85.013 

-5.270 X 10“^ 

unstable 

40.028 




0 

4 


98.416 

-4.276 X 10“^ 

unstable 

40.389 




0 

3 


107.574 

-1.974 X 10-5 

unstable 

40.743 




0 

2 

129.099 

129.406 

-8.392 X 10-^ 

unstable 

41.003 


-0.24 

h 

0 

1 

148.775 

148.579 

-4.386 X 10-’ 

unstable 

42.055 


+0.13 

/4 

0 

0 


165.145 

-2.578 X 10-’ 

unstable 

42.063 




1 

7 


74.783 

+ 1.932 X 10-5 

stable 

39.894 

0.0071 



1 

6 


84.243 

-5.145 X 10-5 

unstable 

39.999 

0.0067 



1 

5 


96.410 

-3.634 X 10-5 

unstable 

40.434 

0.0125 



1 

4 


105.468 

-2.941 X 10-5 

unstable 

40.572 

0.0120 



1 

3 

128.077 

128.176 

-8.586 X lO-'’ 

unstable 

40.999 

0.0141 

-0.08 

/l 

1 

2 

143.649 

143.327 

-8.108 X 10-’ 

unstable 

41.846 

0.0280 

+0.22 

h 

1 

1 


164.511 

-2.909 X 10-’ 

unstable 

42.018 

0.0176 



2 

7 


73.101 

+4.804 X 10-5 

stable 

39.746 

0.0083 



2 

6 

[82.261] 

83.016 

-4.504 X 10-5 

unstable 

39.981 

0.0095 

[-0.92] 

[/lo] 

2 

5 


92.659 

-3.111 X 10-5 

unstable 

40.445 

0.0223 



2 

4 


103.257 

-4.166 X 10-5 

unstable 

40.420 

0.0145 



2 

3 

122.408 

124.141 

-7.354 X 10-® 

unstable 

41.081 

0.0413 

-1.42 

h 

2 

2 

135.160 

135.309 

-3.273 X 10-^ 

unstable 

41.331 

0.0462 

-0.11 

h 

2 

1 


163.248 

-3.474 X 10-’ 

unstable 

41.953 

0.0237 



2 

0 


196.428 

-3.626 X 10-“ 

unstable 

44.808 

0.4159 



3 

7 


71.517 

+9.095 X 10-5 

stable 

39.612 

0.0713 



3 

6 


80.835 

-2.714 X 10-5 

unstable 

40.013 

0.0218 



3 

5 


88.223 

-3.861 X 10-5 

unstable 

40.262 

0.0310 



3 

4 


101.203 

-4.846 X 10-5 

unstable 

40.348 

0.0209 



3 

3 


115.569 

-7.059 X 10-^ 

unstable 

41.145 

0.0822 



3 

2 

133.516 

130.982 

-6.919 X 10-® 

unstable 

41.059 

0.0317 

+ 1.90 

h 

3 

1 

158.789 

160.146 

-4.134 X 10-’ 

unstable 

41.910 

0.0678 

-0.85 

h 

3 

0 

172.214 

171.664 

-2.565 X 10-* 

unstable 

42.888 

0.1774 

+0.32 



degree I and the radial order k of the modes being observed. 
The mode identification appears instead naturally - and some¬ 
what objectively - as a solution (or prediction) of the double¬ 
optimization procedure. Details on the derived mode identifi¬ 
cation and period fit for PG 1219+534 are given in Table 4 (a 
graphical representation of it is also shown in Figure 11). In 
addition to the quantities that reflect the properties of the non- 
radial modes computed from the best-fit model (previously de¬ 
scribed in Subsection 4.2), Table 4 provides the derived distri¬ 
bution of the observed periods (Fobs) as these were matched to 
the theoretical periods (P*)- Again, this distribution is the one 
that minimizes the sum of the squared difference between 
the observed periods and their assigned theoretical periods for 
that model. The relative difference in period AP/P (in %) for 
each pair (Fobs, P*) is also given in this table. 

We find that the simultaneous fit of the nine periods of 
PG 1219+534 obtained with the optimal model retained is ex¬ 
cellent by current standards. The average relative dispersion 
between the fitted periods is AP/P ~ 0.6% and the worst dif¬ 
ference is less than ~ 2%. On an absolute scale, the average 
dispersion between the periods is less than ~ 0.8 s. This accu¬ 


rate and simultaneous fit of nine periods constitutes, by itself, 
a remarkable result in the field of asteroseismology. It is com¬ 
parable in quality (and even slightly better due, in part, to the 
improved tools used to isolate the optimal model) to the results 
achieved by Brassard et al. (2001) for the star PG 0014+067. 
Yet, we emphasize the fact that the accuracy at which the ob¬ 
served periods are measured is approximately one order of 
magnitude better than the mean period dispersion achieved by 
our optimal fit. This suggests that, as good as this fit appears to 
be by current standards of asteroseismology, our equilibrium 
models describing the structure of sdB stars still suffer from 
imperfections that leave significant room for improvement, on 
an absolute scale, in reproducing the current observations. This 
is, of course, one of the goals to pursue in future asteroseis- 
mic studies of sdB pulsators. Interestingly, the marginal detec¬ 
tion of an independent period at 82.261 s (/lo, shown within 
brackets in Table 4), although it was not used at all to con¬ 
strain the search for the optimal model, can indeed, after the 
fact, be associated with a theoretical mode without degrading 
the overall quality of the period fit. This result indicates that 
this period may actually correspond to a real mode excited at 
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very low amplitude in PG 1219+534. In addition, we find that 
the 9 observed periods (10, if we consider /m as real) all fall 
within the predicted band of instability, as it is clearly appar¬ 
ent in Table 4. Hence, the model solution derived from our 
global optimization procedure is also remarkably consistent 
with the prediction of nonadiabatic pulsation theory applied to 
our second-generation models of pulsating hot B subdwarfs. 
The observed periods are identihed with radial {£ - 0) and 
nonradial - 1,2, and 3) /- and p-modes with radial orders 
between k = 0 - 3 (note that the 82.261 s period, if real, would 
correspond to a mode with k - 6, i.e., at the short-period end 
of the band of instability). Interestingly, the mode identification 
inferred from the optimal model presented in Table 4 remains 
unchanged along the line of degeneracy identified in the previ¬ 
ous subsection. Thus, the models found along this line of de¬ 
generacy can indeed be considered as members of a family of 
solutions in this respect. Also of interest, we hnd that the mode 
identification associated with the second (rejected) family of 
solutions found near logg ~ 5.71 is similar to the identihca- 
tion proposed, except that all modes are shifted by Ak = +1. 
This suggests a periodic behavior, with a third family of ac¬ 
ceptable best-fit solutions with modes now shifted by Ak = +2. 
However, such a solution would correspond to a model with a 
surface gravity lower than our bound of logg = 5.65 in our 
search volume in the parameter space, a value in clear conflict 
with available spectroscopic estimates. 

Obviously, from Table 4 and Figure 11, modes of degree up 
to £ - 3 are required in order to account for the mode density 
observed in PG 1219+534. Moreover, all the observed modes 
tend to cluster near the low radial order (low-k) boundary of 
the predicted band of instability. Similar mode distributions are 
actually seen in PG 0014+067 (see Brassard et al. 2001) and 
in a few other well studied EC14026 pulsators. Often, how¬ 
ever, a few periods may be found which correspond to modes 
positioned near the high-k boundary of the instability band, 
thus leaving holes of predicted but unseen modes in the pe¬ 
riod spectrum. A suggestion is that, within the band of excited 
periods, energy powering up the pulsations may be preferen¬ 
tially distributed among the modes of low radial order. These 
modes would then reach observable amplitudes more easily, 
and therefore would be preferentially seen in the light curves 
of EC 14026 pulsators. We also note that we hnd no strongly 
established hierarchy of the mode amplitude as a function of 
£ for the nine modes identihed in Table 4. One would naively 
expect, at least statistically, some trend of amplitude attenua¬ 
tion when £ increases, due to geometric cancellation effects. 
In this context, we note that three of the four dominant pe¬ 
riods are indeed associated with modes of degree £ - 0 and 
1, while periods of lower amplitudes seem to be preferentially 
associated with modes of degree £ - 2 and 3 (with the notice¬ 
able exception of the mode /a identihed as a f = 3 mode). All 
these issues related to mode amplitudes are, however, beyond 
the realm of linear theory and only a nonlinear approach to the 
pulsation phenomenon will be able to address properly these 
questions. In the meantime, a purely observational approach, 
such as searching for the “missing” modes at still higher sen¬ 
sitivity than we could reach with the CEHT, for example using 
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Figure 11. Comparison of the observed period spectrum of 
PG 1219+534 {thick dashed-dotted-line red segments) with the theo¬ 
retical pulsation spectrum of the optimal model. For the latter spec¬ 
trum, solid-line segments indicate excited modes, while dotted-line 
segments correspond to damped modes. All pulsation modes with 
f = 0, 1, 2, and 3 in the period interval 70 - 200 s are illustrated. 
The values of the radial order index k are also indicated for each 
mode. These modes are all acoustic waves (including the /-modes). 
The g-modes have periods that fall outside the range of interest for 
PG 1219+534. 


one nigth of Gemini time, could go a long way toward solving 
this puzzle. 

Einally, we stress that the (k,/) values attributed to the 
periods observed in PG 1219+534 constitute a prediction of 
our seismic analysis which is amenable to independent ob¬ 
servational tests. A promising avenue for independent mode 
identihcation resides in multicolor photometry (as opposed 
to one color, “white light” photometry used, for instance, in 
this study). It is well known from stellar pulsation theory 
that the apparent amplitude of a nonradial oscillation mode 
is a function of wavelength. Moreover, how this amplitude 
changes with wavelength depends on the degree £ of the mode. 
Consequently, measuring relative amplitude ratios of modes 
observed at different wavelength allows, in principle, the iden¬ 
tification of their degree £. Such independent evaluations of 
the geometry of the modes would obviously provide important 
tests of our seismic analyses. Some very encouraging results 
based on that technique have been presented recently by Jeffery 
et al. (2004), including the suggested identification of a mode 
with / = 4 in the hot pulsating sdB star KPD 2109+4401. 
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Table 5. Inferred Properties of PG 1219+534 (V = 13.24 ± 0.20). 


Quantity 


Asteroseismology 


Spectroscopy 


log^ 

Alf (K) 


MJMa 

log(Me„v/M.) 

R/Re(M.,g) 
L/Lo(T,s, R) 

Mv(g, M.) 
d(V, Mv) (pc) 
f'rot (day) 

V,^(R, P„) (km.s->) 
Vgq sin i (km.s“^) 


5.8071 ±0.0057 (0.10%) 
33,640 ± 1,360 (4.04%) 
0.457 ±0.012 (2.63%) 
-4.254 ±0.147 (3.46%) 

0.1397 ±0.0028 (2.00%) 
22.12 ±4.46 (20.2%) 
4.61 ±0.12(2.60%) 
532 ± 37 (6.95%) 
S3.4 
:S 2.1 
2.1 


5.810 ±0.046 (0.79%) 
33,600 ±370 (1.10%) 


22.01 ± 1.85 (8.41%) 
4.62 ±0.06 (1.30%) 
531 ±23 (4.33%) 


$ 10 


4 . 5 . Structural Parameters of PG 1219+534 

The model that has been isolated and which best represents 
the observed properties of PG 1219+534, from both the spec¬ 
troscopy and asteroseismology standpoints, leads to the deter¬ 
mination of the fundamental parameters that define the struc¬ 
ture of this hot pulsating B subdwarf. A first set of primary 
quantities corresponding to the main model parameters is natu¬ 
rally derived from this exercise. These primary parameters are 
the effective temperature the logarithm of the surface grav¬ 
ity logg, the total mass M,, and the mass of the hydrogen- 
rich envelope through the quantity log^(H) ^ log(Menv/4T,). 
A set of secondary parameters then follows from the values ob¬ 
tained for the primary quantities. These secondary parameters 
are the radius R (as a function of M* and g), the luminosity 
L (as a function of Teff and R), the absolute magnitude My 
(as a function of g, and M, in conjunction with the use 
of detailed model atmospheres), the distance from Earth d (as 
a function of V and My). Limits on the rotation period F^ot, 
equatorial rotation velocity yeq(as a function of R and Prot), 
and projected equatorial rotation velocity Veq sin i could also be 
set due to the absence of fine structure in the frequency spec¬ 
trum of PG 1219+534, as multiplets may just be unresolved 
with the time baseline provided by our observations. The val¬ 
ues and limits derived for all these quantities are summarized 
in Table 5. 

An essential ingredient of our asteroseismic study is an 
evaluation of the uncertainties associated with the parameters 

derived for the optimal model. The determination of these un- 

_2 

certainties is closely linked to the behavior of the ^ hyper¬ 
surface in the vicinity of the best-fit solution (see Figures 7 
and 8). Following Brassard et al. (2001), we can model the 
4-dimensional x -function in terms of a quadratic expansion 
considering that the first-order derivatives vanishe at the min¬ 
imum and that the cross-terms in the second-order derivatives 
also vanish, in a first approximation, if we consider the four pa¬ 
rameters involved to be independent. Clearly, that latter state¬ 
ment is not exact, especially in light of the nature of the degen¬ 
eracies that have been uncovered in our search for the best-fit 
solution. These complications will be further discussed below, 
but we will show that they do not, however, significantly affect 
the evaluation of the uncertainties associated with the parame¬ 


ters. This simplified approach leads to the following relation- 

_2 

ship between the variation t+x ^nd the uncertainties, Aa,, on 
the parameters a, (where ai - logg, az = Teft, = M^, and 
04 = log^[H]) 


1=1 


idY 


2 da} 


(Aoi)^ . 


( 2 ) 


The second-order derivatives are evaluated at the position of the 
optimal model and are estimated on the basis of the j^^-values 
found in the grids specifically constructed around this optimal 
model (again, shown in Figures 7 and 8). In order to obtain a 
conservative estimate of the uncertainty on a given parameter, 
we assume that the variation Ax^ is due exclusively to the vari¬ 
ation of that one parameter, thus leading to the relationship 



_2 

The problem now resides in evaluating the quantity Ax ■ The 
merit function that we have defined to measure the quality of 
the period fit is not a standardfor which usual formulae can 
apply to estimate the normal errors. A renormalization of this 
function is needed and can be achieved following the prescrip¬ 
tion of Press et al. (1985; see also Bevington 1969) assuming 
that we have obtained a perfect fit of the observations. In that 
case, the standard x},d at the minimum is equal to the num¬ 
ber of degrees of freedom v, here equal to 5 (9 periods minus 
four free parameters). This leads to a scale factor S between 
the standard minimum x^td ^tid the value of the merit func- 
tion X at the optimal point in parameter space given by 5 = 
5/0.0784 = 63.7755. We next compute the value of that 
must be added to x^i^ to cover a range of parameter space suf¬ 
ficient to reach a given confidence level in the estimates of the 
uncertainties on the parameters. We adopt the Icr limit (68.3% 
confidence level) and compute using the Gammq routine 
of Press et al. (1985) for the truncated gamma function. For the 
case V = 5, we find that = 5.8907. Applying the renor¬ 
malization according to the scale factor S , this corresponds to 
Ax — 0.0924, the value we now use in equation (3). We finally 
find that, at the 68.3% confidence level, the error estimates as¬ 
sociated with the four primary parameters of PG 1219+534 
(i.e., logg, Teff, M„ and log[Menv/M,] ^ log^[H]) correspond 
to the values given in Table 5. These error estimates are also 
illustrated as dotted-line rectangles on the maps provided in 
Figures 7 and 8. Uncertainties on the secondary quantities are 
then easily derived from the values obtained for the primary 
parameters. 

As mentioned in the previous paragraph, this approach re¬ 
lies in part on the assumption that the parameters are indepen¬ 
dent, an assumption that is not well verified in practice due, 
in particular, to the degeneracy in parameter space of the op¬ 
timal model solution. To deal with such a difficulty, one could 
consider, as an alternative method to provide estimates of the 
uncertainties, the relatively strong constraint brought by spec¬ 
troscopy on the determination of the effective temperature. 
Using projection maps of the merit function such as those il¬ 
lustrated in Figures 9 and 10, one can delimit the domain of 
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acceptable values for all model parameters according to the 
limits (i.e., the evaluated errors) imposed by spectroscopy on 
the Teff values. A close examination of these maps indicates, 
however, that the uncertainties derived from this approach are 
of the same order as the values computed from the method in¬ 
troduced by Brassard et al. (2001). We therefore maintain our 
first evaluation of the errors given in Table 5 as representative 
of the Icr uncertainties associated with the derived fundamental 
parameters ofPG 1219+534. 

These asteroseismic results deserve additional comments 
as they are representative of the current potential and limi¬ 
tations of asteroseismology applied to EC14026 pulsators. A 
first significant contribution of asteroseismology to the study 
of sdB stars is the ability to derive values for the surface grav¬ 
ity logg with unprecedented accuracy, i.e., improved by a fac¬ 
tor of ~ 10 compared to current spectroscopic measurements, 
typically. This reflects the fact that the p-mode periods are par¬ 
ticularly sensitive to the gas density which strongly affects the 
speed at which acoustic waves propagate in the stellar interior 
(i.e, the sound speed). And how dense and compact a sdB star 
is largely depends on the log g parameter (see, e.g., Charpinet 
et al. 2002). On the other hand, we find that the asteroseis¬ 
mic value derived for the effective temperature is not well con¬ 
strained, due in part to the fact that the p-mode periods are 
relatively insensitive to this parameter in sdB stars, but also be¬ 
cause solution degeneracies may prevent from uniquely mea¬ 
suring that quantity through asteroseismology only. Notably, 
spectroscopy turns out to be a much more accurate method to 
measure T^g and allows to lift eventual degeneracies in the as¬ 
teroseismic solutions. Therefore, we stress that both methods 
have indeed complementary roles to play in this respect. In ad¬ 
dition to those “classical” parameters, asteroseismology allows 
for the determinations at an interesting level of accuracy of fun¬ 
damental structural quantities, such as the total mass and the 
thickness of the outer H-rich envelope in hot B subdwarfs, two 
key parameters for testing theories of formation and evolution 
of stars on the Extreme Horizontal Branch. These two quanti¬ 
ties cannot generally be determined otherwise, except for some 
masses of sdB stars known through the study of binary systems 
containing an sdB component, but with a relatively poor accu¬ 
racy on the derived value. Eor its part, the determination of the 
mass of the outer H-rich envelope is a pure product of aster¬ 
oseismology and, beside theory, this quantity cannot be con¬ 
strained by other known means. Determinations of MemIMt 
for other EC14026 pulsators through our seismic approach may 
provide, in the future, interesting insight into the internal struc¬ 
ture and evolution of these old extreme horizontal branch stars. 

To conclude this subsection, we point out that we are fully 
aware that our asteroseismic determinations of the global pa¬ 
rameters of PG 1219+534 are only as good, in an absolute 
sense, as the constitutive physics that went into the construc¬ 
tion of the equilibrium models. Euture improvements in our 
ability to understand and model sdB stars at the level of the mi¬ 
crophysics (i.e., the equation of state, the opacity, the radiative 
levitation, initial conditions, and so on) will necessarily lead 
to revised estimates of these parameters and, hopefully, to im¬ 
proved matches of the observed periods. Nevertheless, we do 
point out that solid credibility must already be given to current 


available constitutive physics since, after successfully repro¬ 
ducing the 13 observed pulsation periods of PG 0014+067 to 
better than ~ 0.8% (Brassard et al. 2001), the current models 
are also able to explain the presence of 9 pulsation periods in 
the EC 14026 star PG 1219+534 with an average accuracy of 
~ 0.6% on the periods. 

5. Summary and Conclusion 

We observed the luminosity variations of the relatively bright 
(V - 13.24 ± 0.20), rapidly pulsating hot B subdwarf PG 
1219+534 at high sensitivity at the CEHT during a dedicated 
run in 2004, March. Such observations were part of an ongo¬ 
ing, long-term project to monitor the pulsations of EC14026 
stars at sufficiently high S/N ratios to allow for the detection of 
low-amplitude modes. The objective is to increase the number 
of detected modes in known EC 14026 stars, as it constitutes 
a necessary step for ultimately applying asteroseismic meth¬ 
ods to probe the inner structure of these objects. Our past ob¬ 
serving experience have shown that the “marriage” between the 
Montreal 3-channel photometer Lapoune and the CEHT is par¬ 
ticularly efficient for this purpose. Eor several sdB pulsators 
observed so far in the course of this program, we could typi¬ 
cally double the number of mode detections compared to exis¬ 
tent data from other sites (see, Charpinet 2001). Our most re¬ 
cent observations along this line made no exception to the rule 
and led to the clear identification of nine independent oscilla¬ 
tion periods in PG 1219+534. This constitutes a significant im¬ 
provement over the original data available for this star, which 
allowed for the detection of only four dominant periods (Koen 
et al. 1999). 

Prom the time series analysis using standard methods that 
combine Pourier analysis, least squares fits of the light curve, 
and prewhitening methods, we found that PG 1219+534 has a 
relatively simple pulsation spectrum characterized by well sep¬ 
arated components easily recognizable as independent modes. 
Moreover, we found no indication of close frequency struc¬ 
tures (such as multiplets), at least at the frequency resolution 
achieved during these observations (~ 3.44 //Hz). This indi¬ 
cates that PG 1219+534 is likely a slow rotator, with a lower 
limit on its rotation period that can be set to ~ 3.4 days (i.e., the 
time baseline of our observing run). The contrast is sharp with 
two other members of the EC14026 class, PG 1605+072 (Koen 
et al. 1998) and KPD 1930+2752 (Billeres et al. 2000), which 
are rather exceptional in that both show very complex oscilla¬ 
tion spectra due most certainly to rapid rotation. PG 1219+534 
appears instead as a well behaved sdB pulsator with a rich - yet 
simple enough - period spectrum ideally suited for detailed as¬ 
teroseismology. Interestingly, we found that, for this particular 
star, one single night of excellent CEHT data turned out to have 
all the seismic information necessary for a subsequent success¬ 
ful asteroseismological exercise. Our detailed analysis of the 
light curves of PG 1219+534 also illustrated the fact that, for 
EC 14026 pulsators having an easily resolved period spectrum, 
the search for a better signal-to-noise ratio is of higher interest 
for the sake of asteroseismology than a longer time baseline 
and/or coverage. Indeed, we showed that mixing high quality 
light curves with low S/N ratio data (due to degraded observ- 
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ing conditions) can basically ruin all the benefit of having the 
high sensitivity data, since the overall noise level in the Fourier 
domain will be dominated by the noise introduced by the worst 
time series included. This has, of course, implications on how 
one should treat heterogeneous data sets such as those gathered 
during multisite campaigns, for instance. We mention, along 
this line, that the best asteroseismic return for more complex 
EC 14026 pulsators that need multisite contiguous coverage to 
unambiguously resolve their pulsation spectrum would require 
the acquisition of high sensitivity photometric data which is as 
homogeneous as possible in quality. In this context, we are en¬ 
gaged in the organisations of bi- (or tri-) site campaigns that 
would involve 4m-class telescopes positioned at strategic lon¬ 
gitudes. Our efforts are presently directed toward a joint utili¬ 
sation of the CFHT with Lapoune and the WHT with Ultracam 
as both telescopes are ideally located geographically and would 
provide an almost complete coverage. This, of course, would 
be in addition to recent efforts to perform multicolor fast pho¬ 
tometry (Jeffery et al. 2004) and time resolved spectroscopy 
(O’Toole et al. 2004) on these pulsating stars. 

Based on the nine oscillation periods identified in the light 
curve of PG 1219+534, we have attempted a detailed aster¬ 
oseismic analysis of this EC 14026 star. Our approach relied 
on the well-known forward method with the goal of finding 
objectively the model that would best match the set of pe¬ 
riods observed in this star with a set of theoretical periods. 
We have used the second generation sdB models of Charpinet 
et al. (1997) to compute the theoretical periods since those have 
proved to account quite well for the class properties of the 
EC14026 pulsators (Charpinet et al. 2001). Due to the lack of 
mode identification from the available “white light’’ photomet¬ 
ric data, we performed a double-optimization procedure that 
takes place simultaneously at the period matching level and in 
the model parameter space, leading objectively to the best fit 
of the asteroseismic observations. In the present study, we have 
used the Toulouse package of numerical tools specifically de¬ 
veloped, over the years, for the purpose of asteroseismology 
and which includes a genetic algorithm based period matching 
code (dedicated to the first optimization step), a parallel ge¬ 
netic algorithm based code, and a parallel grid computing code 
allowing for an efficient exploration of the vast parameter space 
(the second optimization step). 

Our thorough exploration of parameter space led us to first 
isolate two families of model solutions that equally best-match 
the observed periods of PG 1219+534. These families of mod¬ 
els have surface gravities around logg ~ 5.81 and logg ~ 5.71, 
respectively, and correspond to slightly different mode identifi¬ 
cations. This degeneracy in the asteroseismic solution could not 
be lifted solely on the basis of period fit quality, and additional 
constraints were required to select the “correct” one. Such con¬ 
straints were provided by the available spectroscopic measure¬ 
ments for the atmospheric parameters of PG 1219+534. An 
ongoing project aimed at improving the spectroscopic char¬ 
acterization of sdB stars (including both the rapid and slow 
pulsators), and which is based on the acquisition of high res¬ 
olution, high S/N ratio spectra with the blue spectrograph at 
the new MMT, has proved essential in this context. The inde¬ 
pendently derived atmospheric parameters for PG 1219+534, 


Teff = 33,600 + 370 K and logg = 5.810 + 0.046, allowed 
us to clearly reject the family of model solutions found at low 
surface gravity (i.e., near logg ~ 5.71). Concentrating on the 
preferred solution located near logg ~ 5.81, we found however 
that another kind of degeneracy exists in the solution which, 
again, precludes from uniquely isolate an optimal model from 
asteroseismology alone. This problem occurs due to concur¬ 
rent effects of model parameters on the pulsation periods, lead¬ 
ing, in the worst case, to a line-degeneracy of the best-fit so¬ 
lutions that mostly link the total mass parameter with the ef¬ 
fective temperature parameter. We found, indeed, that seismic 
solutions of comparable quality exist for all values of the total 
mass of the star, but these solutions have values for that 
are correlated. The trend uncovered is monotonic and indicates 
that optimal models of low-mass (down to 0.30 Mq) require 
cooler effective temperatures, while optimal models of high- 
mass (up to 0.53 Mq, but also verified with a few additional 
calculations at higher masses) need hotter effective tempera¬ 
tures. It is again spectroscopy that provided the additional in¬ 
formation needed to bypass this difficulty, as using the well- 
constrained value of Teff derived from this method led to the 
determination of a unique viable optimal model solution that 
fit consistently PG 1219+534, both asteroseismologically and 
spectroscopically. 

The basic properties of PG 1219+534 as inferred from our 
combined spectroscopic and asteroseismic approach have been 
summarized in Table 5. The quest for an optimal seismic model 
of PG 1219+534 has revealed that accurate spectroscopy re¬ 
mains an essential ingredient for detailed asteroseismology of 
EC 14026 stars in order to converge toward a unique solution. 
A relatively precise determination of the surface gravity us¬ 
ing spectroscopic techniques is generally helpful to discrimi¬ 
nate between families of potential seismic solutions. Once the 
“right” family has been chosen, though, the constraint brought 
by asteroseismology on the value of log g is much tighter than 
the precision that spectroscopy can currently provide. This is 
one of the outstanding outcomes of asteroseismology as it leads 
to measurements of the surface gravity of EC 14026 stars at 
unprecedented accuracy (a gain of a factor ~ 10, typically). 
The opposite is true for the effective temperature and the spec¬ 
troscopic value for this parameter must clearly be preferred. 
Moreover, in the case of PG 1219+534, the line degeneracy 
uncovered for the seismic solution precludes an independent 
measurement of Teff through asteroseismological means. 

Another outstanding outcome of asteroseismology is the 
measurement of the total mass and the mass of the H-rich 
envelope at an interesting level of accuracy. In our study of 
PG 1219+534, we found however that the total mass could not 
be inferred independently of the spectroscopic measurement 
of the effective temperature. Nonetheless, the value derived, 
M = 0.457 ± 0.012 Mq, indicates a mass of PG 1219+534 
that is close to the canonical value for sdB stars. Although this 
value of the mass depends on how reliable the spectroscopic 
estimate of is (particularly in regards to eventual system¬ 
atic effects), we stress that significantly different masses than 
the derived - close to canonical - value, such as some scenar¬ 
ios of binary evolution suggest, would require a shift of sev¬ 
eral thousand Kelvins in effective temperature to be compatible 
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with asteroseismology. Such a drastic change in the evaluation 
of Teff from spectroscopy is highly unlikely and, therefore, we 
are confident in concluding that PG 1219+534 has a mass close 
to the canonical mass of extreme horizontal branch stars. 

We found that our optimal model solution is able to repro¬ 
duce simultaneously the 9 periods observed in PG 1219+534 
with an average dispersion of only ~ 0.6%. The observed peri¬ 
ods correspond to low-order { - Q,\, 2, and 3 acoustic modes 
that are, indeed, predicted to be excited according to nona- 
diabatic pulsation theory. It is a remarkable result (not often 
achieved so far in the field of asteroseismology) that the model 
solution uncovered for PG 1219+534 can, at the same time, re¬ 
produce all the periods observed in this star at a relatively high 
level of accuracy, be consistent with nonadiabatic theory, and 
satisfy the spectroscopic constraints. Such a consistency be¬ 
tween three independent aspects of the modeling of these pul¬ 
sating stars is a result that was not guaranteed at the outset. This 
is strong indication that the basic constitutive physics which 
enters in the construction of our current models used for the as- 
teroseismic analyses of EC 14026 stars is sound, especially as 
it follows similar results already obtained for another rapid sdB 
pulsator, PG 0014+067 (Brassard et al. 2001). Furthermore, 
the excellent agreement between observations and theory at the 
nonadiabatic level brings further proof that the iron bump opac¬ 
ity mechanism of Charpinet et al. (1997) is indeed at the origin 
of the EC 14026 phenomenon. 

Finally, we note that the mode identification resulting from 
our detailed asteroseismic analysis of PG 1219+534 consti¬ 
tutes a prediction that, in principle, can be independently tested 
using multicolor photometry. Along this line, we initiated a 
project to use the Far UV Spectroscopic Explorer (FUSE) 
as a fast-photometer (a possibility offered by the so-called 
“Time-Tag” mode of the instrument) that would provide a light 
curve of PG 1219+534 in the FUV (near ~ lOOOA) bandpass. 
Combined with nearly simultaneous ground based observations 
in the optical, this project should lead to strong independent 
constraints on the ( value of the modes seen in PG 1219+534, 
from which the predicted mode identification reported in this 
paper will be checked. 
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